NASA Technical Memorandum 83253 

NASA-TM-83253 19840020449 


AN IMPLEMENTATION OF THE DISTRIBUTED PROGRAMING 
STRUCTURAL SYNTHESIS SYSTEM (PROSSS) 


JAMES L. ROGERS^ JR. 


DEC 2 4 1981 


f:Al-1bLL> 

LIBRARY, NASA 
HAWPT9N, vmoiM'A 

DECEMBER 1981 


FOR I 


JLY DOMESTIC DISSEMIN^ON 


Because oif its significant early^Tnimercial potential, this 
informationWhich has been c^eloped under a U.S. Gov- 
ernment proOTam, is beino4^seminated within the United 
States in advar^ of gederal publication. This information 
may be duplicaf^Wa^^^ed by the recipient with the ex- 
press limitation not be published. Release of this 

information to o^lra^estic parties by the recipient shall 
be made subjj^To th^j^limitations. 

Foreign re^ase may bel^^de only with prior NASA ap- 
proval a^ appropriate ex^^ licenses. This legend shall 
be m^^ed on any reproductio^f this information in whole 
or jj^art. 

nber 31 , 1983 


aview for general release 


IWNSA 

National Aeronautics and 
Space Administration 

Langley Research Center 

Hampton, Virginia 23665 



AN IMPLEMENTATION OF THE DISTRIBUTED PROGRAMING 
STRUCTURAL SYNTHESIS SYSTEM (PROSSS) 

James L. Rogers, Jr. 

NASA Langley Research Center 
Hampton, Virginia 
INTRODUCTION 


Numerous approaches have been documented for combining optimization 

I 

techniques ^ith an analysis capability (e.g., refs. 1-4). The PROgraming 

I 

structural Sy>^thesis System (PROSSS), described in reference 4, was first 

implemented entirely on a CDC CYBER computer. This implementation is 

described in reference 5. PROSSS was then distributed between a mainframe and 

a minicomputer (ref. 6) to provide the user easier access to intermediate 

results, more control over the flow of the problem, and an interactive 

modeling and data generation capability (ref. 7). 

This document describes the implementation of the distributed version of 
PROSSS, the second step in a series of implementations. After an overview 
which explains PROSSS in general with respect to this implementation, the 
remainder of this document is devoted to describing (1) each component of the 
system, (2) how to examine the intermediate results, (3) how to restart the 
system, and (4) results from a sample problem exercising the options of this 
implementation. Appendices contain listings to clarify the descriptions' of 
the components of the system contained in the body of this paper. 

It is not the purpose of this paper to be a self-contained user's guide 
for PROSSS. Its purpose is to demonstrate a method for implementing a 
flexible software system that combines large complex programs with small. 



user-supplied, problem-dependent programs and distributes their execution 
between a mainframe and a minicomputer. 

OVERVIEW OF THIS IMPLEMENTATION 

This distributed implementation of the PRO graming _Structural Synthesis 
System (PROSSS) combines a general purpose finite-element computer program for 
structural analysis (SPAR, ref. 8), a state-of-the-art optimization program 
(CONMIN, ref. 9), and several user-supplied, problem-dependent computer 
programs. The results are flexibility of the optimization procedure 
organization and versatility of the formulation of constraints and design 
variables. 

Figure 1 is a flow chart of the analysis-optimization process for this 
implementation. The process results in a minimized objective function, 
typically the mass. Notice that the analysis and optimization programs are 
executed repeatedly by looping through the system until the process is stopped 
by a user-defined termination criterion. This is referred to as the 
repeatable part of PROSSS. However, some of the analysis, such as model 
definition, need only be done one time and the results are saved for future 
use. This analysis is performed outside of the loop and referred to as the 
nonrepeatable part of PROSSS. The user must write some small, simple FORTRAN 
programs to interface between the analysis and optimization programs. One of 
these programs, the front processor, converts the design variables output from 
the optimizer into a suitable format for input into the analyzer. Another, 
the end processor, retrieves the behavior variables and, optionally, their 
gradients from the analysis program and evaluates the objective function and 
constraints and optionally their gradients. These quantities are output in a 
format suitable for input into the optimizer. These user-supplied programs 
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are problem-dependent because they depend primarily upon which finite elements 

. are being used in the model. 

The analysis is always performed on the mainframe, in this implementation 

If- 

a CDC 6600, because the CPU on the mainframe is much faster than that of the 
minicomputer. The end processor is also executed on the mainframe to reduce 
the large amounts of data generated by the analysis program to the small 
amount of data (e.g., objective function and constraints) required by the 
optimizer. This limits the amount of data transmitted between the two com- 
puters. 

The optimization is always performed on the minicomputer, in this imple- 
mentation a PRIME 750, because the optimizer requires a significant memory 
allocation and the minicomputer has virtual memory. The front processor is 
also executed on the minicomputer and, as in the case of the end processor, 
the amount of data transmitted between the two computers is reduced to the 

* minimum. The data transmitted from the front processor to the analysis pro- 
gram consist of updated design variable information. 

Five options exist for organization of optimization procedures comprising 
nonlinear and piecewise linear programing with analytical and finite- 
difference gradients. These options are shown in Table I. One of the primary 
differences between the distributed implementation of PROSSS and the first 
implementation described in reference 5 is the method of controlling the flow 
of the options. In the first implementation, each option is controlled by a 
complex sequence of Control Data Corporation-Cyber Control Language (CDC-CCL 
(ref. 10)), while the flow of the options in the distributed implementation of 

i 

PROSSS is controlled by FORTRAN programs using FORTRAN-callable system sub- 

* routines (ref. 11). Although presented in the context of structural analysis. 
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these organizations are independent of the type of analysis. The same system 
concept could be used for aerodynamic optimization of, for example, a wing, if 
the analyzer had a capability for computational aerodynamics. 

The system is intended to be used in the following three basic ways: 

1. as a research tool for development of optimization techniques where 
it offers an interface to efficient analysis and flexibility of 
execution and sequencing, 

2. as a test bed for trying new structural analysis techniques tailored 
toward efficient use in optimization loop, and 

3. as an application tool that can be adopted to a very wide scope of 
different types of problems. 

Before trying to implement this system, the user should be thoroughly 
familiar with references 4, 5, 6, 8 and 9. 

COMPONENTS OF THE SYSTEM 

There are seven primary components for the distributed implementation of 
PROSSS including (1) the driver program, (2) the front processor, (3) CONMIN- 
the optimizer, (4) procedures transmitting jobs between the PRIME minicomputer 
and CDC mainframe computer, (5) SPAR-the analyzer, (6) the end processor, and 
(7) the plotting program. Each component and the files associated with them 
are explained in detail in the following subsections. Files generated by the 
system are not shown. Only those files that must be created by the user are 
shown. 

The Driver Program 

As previously mentioned, each option in PROSSS is controlled by a FORTRAN ^ 

program using FORTRAN-callable procedure files. While each option uses its 
own main driver prograin, there is a group of eight subroutines common to each 
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option. A typical flowchart of an option is shown in figure 2. The listings 
of the program and subroutines are shown in Appendix A. The driver program 
performs the following functions: (1) initialization, (2) restarting, (3) 

controlling the flow of the option, and (4) termination. 

The driver program first reads one line from a file named NAMEij where ij 

is the option number without the decimal. For example, a file named NAMEll 
contains the input data for option 1.1. This first line in NAMEij contains 
the following variables: 

NRSTRT (8 - not restarting 

1 - restarting 

ICR 0 - data not created by a data generator 

1 - data created by a data generator 

NR 0 - no nonrepeatable analysis is to be done 

1 - perform nonrepeatable analysis 

NOPT the option number (e.g. 11) 

NAMES the name of the file to be sent to the CDC 6600 (four 
characters in length) 

COUNT the number of passes made through the analysis program 
(initially one) 

The format is (3(I1,1X) ,I2,lx,A4,lX,I3). 
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An example of a NAMEij file is as follows: 


0 0 0 11 

ANOS 

1 

STRP3 


file containing 

CONPl 


file containing 

CONS 


file containing 

RRSNG 


file containing 
gradients 

JIM 


UFD name 

DPROSS 


SUBUFD name 

SSPROUT 


file containing 

SCNMNIO 


file containing 

NRLANG 


file containing 

123456 


user number 

XYZ 


password 

123456 


charge number 

EPFUBIP 


file containing 

EPFUIN 


file containing 

FPFUS 

0 

file containing 

CONMSl 

1 

file containing 
compile) 


starting design variables 
CONMIN parameters 
constants for front processors 
repeatable runstream with no analytical 


SPAR output 

end processor output 

non-repeatable SPAR library 


end processor 

input to end processor 

front processor (0 mean compile) 

CONMIN main program (1 means do not 


If NRSTRT is equal to , zero, subroutine INIT is called to create a 
temporary procedure file called INITPRC. INIT reads the remainder of file 
NAMEij and creates edit procedures on file INITPRC to replace general file 


names in procedure files executing other components of the system with 


specific file names depending upon the problem being solved. Files for other 
components of the system are created depending on the option number. (This is 
transparent to the user.) 
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If ICR is not equal to zero, then INIT calls subroutine CRSNPT. This 
subroutine reads data from a file called NAMECRRS, edits a file created by the 
data generator, and inserts a file containing problem-dependent data. An 
example of a NAMECRRS file is as follows: 

1 number of design variables 

SPARIN file created by a data generator for input to SPAR 

TOPSPARIN file containing problem dependent data for input to SPAR 

If NR is not equal to zero, then INIT calls subroutine NONRPT. This sub- 
routine reads data from a file called NAMENR and creates a procedure file 
(TONOSNR, described below) to be sent to the CYBER to perform the nonrepeat- 
able part of the analysis. An example of a NAMENR file is as follows: 

TSPARIN non-repeatable SPAR runstream 

700Ji!)0 field length required on CDC 6600 

NONl non-repeatable job name 

If NRSTRT equals 1, then subroutine RSTRT is called instead of subroutine 
INIT. This subroutine sets up a procedure file named RSTRTPRC to reinitialize 
the files saved for restarting the program. 

If the nonrepeatable part of the analysis is required (NR = 0), a job has 
been sent to the CDC 6600 to perform this analysis. At this point, subroutine 
STCHEK is called. This subroutine calls subroutine COMSTA which returns the 
status of a file that has been transmitted to the CDC 6600 computer. If 
COMSTA determines that the job has not yet completed executing on the CDC 6600 
(because no file has been returned to the PRIME), then it returns a status of 
zero. If the job has completed executing on the CDC 6600 and a file has been 
returned to the PRIME, the status is returned as one. When the status cannot 
be determined, an error has occurred and the status is returned as -1. 
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Subroutine STCHEK checks the status returned from COMSTA. If the status is 
one, then control is passed back to the driver program. If the status is -1, 
subroutine ERRMSG is called to print an error message, and the job is termi- 
nated. If the status is zero, a system routine is called to put the job on 
the PRIME to "sleep" for 2 minutes while waiting for the CDC 6600 to finish 
processing. After 2 minutes, the job "awakens" and again calls subroutine 
COMSTA to determine if the file has been returned from the CDC 6600, The 
"sleep" and "awake" process continues until the CDC 6600 has finished pro- 
cessing the analysis and returned the output file to the PRIME. This concept 
makes efficient use of the PRIME'S resources. 

Once the nonrepeatable part of the analysis program has been completed, 
the driver enters the optimization loop shown in figure 1. It will remain in 
this loop until an error is encountered, a loop limit (default is 999) is met, 
or the termination criteria are satisfied. In the optimization loop, proce- 
dure files executing the optimization program and front processor program 
(described below) are called. The front processor creates a file containing 
updated design variable information to be transmitted to the CDC 6600. The 
file sent to the CDC 6600 to perform the repeatable part of the analysis must 
have a unique job name to facilitate the return checking done by subroutine 
COMSTA. Subroutine UPDSF is called to create a temporary procedure file 
(UPDTE) for assigning this unique job name. The job name is created from the 
NAMES variable from the first line of the NAMEij file coupled with the loop 
counter. 

Two temporary files containing the names of files pertaining to the 
specific problem at hand are created during initial i zation from two permanent 
files containing general file names. One temporary file, the names coming 
from the NAMES variable on the NAMEij file, is created from a permanent file 
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called TONOS. TONOS is a file containing CDC-CCL statements to execute the 
repeatable analysis (SPAR) and end processor on the CDC 6600. The second tem- 
porary file, called TSENDPRC is created from the permanent file SENDPRC. 
SENDPRC is a PRIME procedure file for appending the data file created by the 
front processor to the temporary file created from TONOS and then transmitting 
that file to the CDC 6600. Once the TSENDPRC file has been executed and 
processing is transferred to the CDC 6600 for the repeatable part of the 
analysis, then subroutine STCHEK is again called to test for the status of the 
file to be returned from the CDC 6600. The same checking, "sleeping," and 
"awakening" pattern described above for the nonrepeatable analysis is also 
followed here. 


4 
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A sample listing of a TONOS file is shown below. This listing contains the 
general names that are replaced by specific, problem dependent names during 
initialization. 


JOB,Tl!3!9,CM700i9!3. 

USER,USERNO,PASSWRD. 

CHARGE, CHARGENO.LRC. 

DELIVER. $UFD>SUBUFD>TCMNIO 
RFL, 70000. 

REDUCE(-) 

MAP, OFF. 

COPYCR, INPUT, SPFPOUT. 

REWIND, SPFPOUT. 

GET,RUNSTREAM,MERGFP. 

EDIT,RUNSTREAM,,MERGFP. 

EDIT,RUNSTREAM,,MERGFP,EDOUT. 

REWIND, RUNSTREAM. 

GET,SPAR=SPAR14I,DCU=DCU14I. 

GET,SPARLA=SPAROUT. 

SPAR, RUNSTREAM, SPAROUT. 

REWIND, SPAROUT. 

REWIND, SPARED. 

GET,EPXXXX,ENDINXX. 

GET,SPARLIB. 

LDSET(LIB=SPARLIB,PRESET=ZERO) 

EPXXXX,ENDINXX,CNMNIO. 

REPLACE, SPAROUT=DUMOUT,CNMNIO=DUMIO. 
REWIND, CNMNIO. 

C0PYSBF,CNMNI0, OUTPUT. 
DAYFILE(L=DISOUT) 

REPLACE, DISOUT. 

EXIT. 

DAYFILE(L=DISOUT) 

REPLACE, DISOUT. 



A sample SENDPRC file is; 

SAVE TEMP 

DELETE UPDTE 

APPEND SENDFILE SPFPOUT 

DELETE CNMNIO 

RJSEND SENDFILE CDC 

DELETE SPFPOUT 

R TEMP 

After the analysis output file, TCMNIO, is returned from the CDC 6600, 
the driver program checks for the existence of a file called GONOGO. GONOGO 
is created by the CONMIN main program when the termination criteria have been 
met. If GONOGO exists, then the optimization process stops. If GONOGO does 
not exist, then the optimization loop continues by calling the EDPRC procedure 
file. This procedure file edits TCMNIO to remove any extraneous information, 
such as the dayfile, returned from the CDC 6600. This leaves only the 
objective function and constraint data required for input into CONMIN. At 
this point, the optimization loop begins a new iteration. A sample EDPRC file 
is: 


SAVE TEMP 
ED TCMNIO 
L END CNMNIO 
01000 
T 

L BEGIN CNMNIO 
D 

UNLOAD CNMNIO 1000 
FIL 

DELETE TCMNIO 
R TEMP 
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The Optimizer 

The procedure file called from the driver program to execute CONMIN, the 

optimizer, is named CNMNPRCl. Whenever any procedure file is called from the 

driver program, the first instruction in the procedure file is to save the 

current address. The final instruction in the procedure file is to return 

control to the saved address. The CNMNPRCl procedure file performs four 

functions: (1) saves restart files, (2) opens and closes files for CONMIN 

input and output, (3) executes CONMIN, and (4) save cumulative plotting and 

optimization data. A listing of CNMNPRCl is below: 

SAVE TEMP 
ED PASS 
FIL SPASS 
ED CNMNIO 
FIL SCNMNIO 
. ED CONREST 
FIL SCONREST 
OPEN TPLOT 12 2 
OPEN PLTDATA 11 3 
OPEN PASS 7 3 
OPEN CNMNIO 5 3 
OPEN CONREST 33 
OPEN CONOUT 2 2 
OPEN CONPAR 1 1 
SE6 CONMIN 
C 1 
C 2 
C 3 
C 5 
C 11 
C 12 

APPEND PLTDATA TPLOT 
DELETE TPLOT 

APPEND CONOUTHOLD CONOUT 
DELETE CONOUT, 

R TEMP, 

r-/ 
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Three files are required to restart CONMIN. File PASS contains the 
number of passes made through C0NMIN--one initially and two on each subsequent 
pass. File CNMNIO contains the objective function and constraints found in 
the repeatable analysis. File CUNREST contains all of the remaining data 
saved from the preceding pass through CONMIN. These files are saved as 
temporary files named SPASS, SCNMNIO, and SCONREST, respectively. 

Other files open for CONMIN input and output include TPLOT, PLTDATA, 
CONOUT, and CONPAR. TPLOT contains the plotting data created in this (and 
only this) pass through CONMIN. PLTDATA contains a cumulative history for the 
objective function, design variables, and constraints for the entire run. 
CONOUT contains the output listing for only this pass through CONMIN. CONPAR 
contains the output listing for this pass through CONMIN. CONPAR contains the 
initialization values for the parameters input to CONMIN. This file must be 
created by the user and is input via a free-field format. A description of 
these variables, except for JSC, can be found in reference 9. The variable 
JSC is used to determine the linear constraint identification vector, JSC. If 
JSC is -999, then all of the constraints are known to be a linear function of 
the design variables. If JSC is zero, then all of the constaints are 
nonlinear. If JSC is other than these two numbers, then a file, ISCFILE, must 
be created by the user containing the constraint numbers (in 14 format) of the 
linear constraints if JSC is positive or the constraint numbers of the 
nonlinear constraints if JSC is negative. JSC in this last case is the number 
of constraints, in ISCFILE. ISCFILE is opened and closed within CONMIN. 

Another file, STARTX, is also opened and closed within CONMIN and must be 
created by the user. This file contains the number of design variables in 15 
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format followed on succeeding lines by the starting values of the design 
variables in 6E13.5 format. A sample CONPAR file is: 

5 3 20 190 0 1 0 

4 0 0 0.0 0.0 0.0 0.025 0.0 0.0 

0.125 5. .05 0.0000000001 0 20 200 100 100 200 

0.0 0.0 0 0.005 0.1 1.0 1.0 10. 10. 

A sample STARTX file is: 

3 

0.50000E=0.1 0.10000E 01 0.40000E 01 

Once the restart files have been saved and the necessary files have been 
opened, the CNMNPRCl executes CONMIN. CONMIN consists of a group of nine 
subroutines as defined in reference 9. These routines are never changed by 
the user. The routines are called by a CONMIN main program written by the 
user dependent upon the problem to be solved. An example of a main program is 
listed in Appendix B. In addition to calling the CONMIN subroutines, the 
CONMIN main program is responsible for reading and writing data from and to 
the above files. After executing CONMIN, CNMNPRCl closes all files that 
remain open. 

The final task performed by CNMNPRCl is to append the temporary files 
containing the plotting data (TPLOT) and CONMIN output listing (CONOUT) to the 
permanent cumulative files PLTDATA and CONOUTHOLD, respectively. TPLOT and 
CONOUT are deleted and control is returned to the driver program. 

The Front Processor 

The procedure file called from the driver program to execute the front 
processor, FPFUS, is named FRONTPRC. The front processor is supplied by the 
user and converts the design variables output from CONMIN on file CNMNIO to 
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a file SPFPOUT formatted suitable for input into the SPAR analysis program. A 

listing' of FRONTPRC follows: 

SAVE TEMP 
OPEN CNMNIO 1 1 
OPEN SPFPOUT 2 2 
OPEN CONSXX 3 1 
R *FPROCESSOR 
C 1 
C 2 

RTO 

The FRONTPRC procedure file opens the CNMNIO and SPFPOUT files. In 
addition, a file called CONSXX, supplied by the user if necessary and 
containing certain constants such as the cross section dimensions of a beam, 
is also opened. The format of the CONSXX file is determined by the user 
because the user writes the front processor program dependent upon what 
elements are used in the finite-element model. An example of a front 
processor is then executed, all open files are closed, and control is returned 
to the driver program. A sample CONSXX file for the beam parameters B10, B20 
and T0 follows: 

12. 30. 2. (beam parameters) 

The Analyzer— Nonrepeatable Part 

The INITPRC procedure file (described above in the section on the driver 
program) transmits the TONOSNR procedure file to the CDC 6600 to execute the 
nonrepeatable part of the analysis. The nonrepeatable SPAR runstream (NRRS) 
is appended to the TONOSNR file for transmittal to the CDC 6600. NRRS is 



copied to a local file and SPAR is then executed using NRRS as input. The 

results are saved on a SPAR library file (NRLA) for later use in the 

repeatable part of the analysis. A file named TNREPT is returned to the PRIME 

indicating that the nonrepeatable analysis has completed execution. A listing 

of TONOSNR follows. 

JOB.TIOO.CMFLX. 

USER.USERNO.PASSWRD. 

CHARGE, CHARGEI\!0,LRC. 

DELIVER. $UFD>SUBUFD>TNREPT 
MAP, OFF. 

COPYCR, INPUT, NRRS. 

REPLACE, NRRS. 

.* THE PROCEDURE CREATES A SPAR LIBRARY 

.* FROM THE NON-REPEATABLE PART 
* 

• 

GET,SPAR=SPAR14I,DCU=DCU14I. 

GET, NRRS. 

RFL,FLX. 

REDUCE,-. 

SPAR, NRRS, NSPROUT. 

REPLACE, SPARLA=NRLA. 

!* TEST TO SEE IF GRADIENTS ARE REQUIRED 

★ 

• 

IFE,(NR0PT.EQ.13.0R.NR0PT.EQ.23) , GRADIENTS. 

GET,INPT=INPTXX. 

GET,EDIT1,EDIT2,BLDELDB. 

REWIND, NRRS. 

.* 

.* EDIT OUT ALL BUT ELD INPUT IN RUNSTREAM 

EDIT, NRRS, EDIT1,ED0UT. 

REWIND, NRRS. 

.* 

.* CREATE SPAR RUNSREAM TO FIND DERIVATIVES 

BLDELDB,INPT,NRRS,RSOUT. 

REWIND, RSOUT. 

EDIT,RS0UT,,EDIT2,ED0UT. 

REWIND, RSOUT,, SPARLA. 

.* EXECUTE SPAR AGAIN TO FIND DERIVATIVES OF 
.* THE STIFFNESS MATRIX WITH RESPECT TO 

.* THE DESIGN VARIABLES 

SPAR, RSOUT, NSPROUT. 

RETURN, INPT,EDIT1,EDIT2,BLDELDB,ED0UT. 

RETURN, RSOUT, TAPE21,TAPE22, NSPROUT. 

REPLACE, SPARLB-NRLA. 

ENDIF, GRADIENTS. 

RETURN, SPARLB=NRLA. 

ENDIF, GRADIENTS. 

RETURN, SPAR, DCU, SPARLA, NRRS. 

DAYFILE(L=DISOUT) 

REPLACE, DISOUT. 

EXIT. 

DAYFILE(L=DISOUT) 

REPLACE, DISOUT. 
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An example of a non-repeatable SPAR runstream (NRRS) is as follows: 

CXQT TAB 

START 80# 5$ ROTATIONS ABOUT Y EXCLUDED 
TITLE” FUSELAGE M0DEL#FSPAR1 

TEXT 

” MEMBRANE-ROD-BEAM FUSELAGE MODEL 
"NONREPEATABLE PART 

JLOC$ FUSELAGE DIA. 800. CM .# LENGTH=800 . CM. 
F0RMAT»2$CYLINDRICAL COORDINATES 

1 AOO. 0. 0. AOO. 337,5 0. 16 1 5 
16 AOO. 0. 800. AOO. 337.5 800. 

MREF 

FORMAT*? 

1 -2 0 . 0 . 10000000 . 

2 1 0 . 0 . 10000000 . 

MATC 

1 .72+6 0.3 .0028 22. -6t AL-ALLOY# METRIC UNITS 

E23 SECTION PROPERTIES $ROD ELEMENTS 

1 A.168$AREA OF THE RODS 

E21 SECTION PROPERTIE.SSBEAM ELEMENTS 

DSY 1 1680A. 0. 1262.7 0. 108. l^A. 0. 6.078A 0. 0. 

0. 0. 0. -8.7778 17. 3.2222 17. 3.2222 -17. -8.7778 -17. 
SHELL SECTION PROPERTIES 
I O.ISSKIN THICKNESS 
CONSTRAINT CASE 1 

ZERO 1#2#3;1#16$ CANTILEVER THE FUSELAGE 


CXQT ELD 


tSTART 

A1 A2$ 

E23$R0D ELEMENTS 

^ A2 A3$ 

NSECT«1$ 

/f A3 AA$ 

NREF=2 

X AA A5$ 

1 17 1 A 3 1 $ 

X A5 A6$ 

A 20$ 

A6 A7$ 

5 21$ 

A7 A8$ 

6 22$ 

A8 33$ 

52 68$ 


53 69$ 

y 

5A 70$ 

/ $END 

7 23 1 A 10 1$ I 

$START 

$EN0 1 

FA1$ MEMBRANE PANELS 

$START \ 

NSECT*1$ 

F21 

1 17 18 2 2 16 1 $ 

NSFCT»1$ ' 

A9 65 66 50 2 16 1$ 

NREF*1 

17 33 3A 18 1 1 2 2 16$ 

1 2 2 16 2 16 $ 

23 39 AO 2A 1 1 9 2 16$ 

A9 50 2 16 2 16$ 

32 A8 33 17$ 

33 3A$ 

A8 6A A9 33$ 

3A 35$ 

$END 

39 A0$ 

CXQT DCU 

AO Al$ 

TOC 1 


CXQT EXIT 


17 



The Analyzer--Repeatable Part 


The SENDPRC procedure file transmits the TONOS procedure file (both files 

discussed in the above section on the driver program) to the CDC 6600 to 

execute the analysis (SPAR) and end processor (EPXXXX) programs. The SPFPOUT 

file created by the front processor is appended to the TONOS file for 

transmittal to the CDC 6600. The SPAR library file (SPARNAME) containing data 

created during the nonrepeatable analysis and the repeatable SPAR runstream 

file (RUNSTREAM) both reside as permanent files in the CDC 6600 disk storage. 

SPFPOUT is copied onto a temporary file on the CDC 6600 and merged into 

RUNSTREAM file using a file named MERGFP consisting of CDC TEXT EDITOR (ref. 

12) commands. SPAR is then executed using the RUNSTREAM and SPARNAME files as 

input. The objective function and stresses are computed in SPAR and saved on 

a SPAR library file (SPARED). This file is then used as input to the end 

processor. A listing of the MERGFP file is as follows: 

MERGE :/SPFP0UT/,/$ MERGE NEW/;1 
R 

RS:/E+/,/+/;* 

RS:/E-/,/-/;* 

END 
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An example of a repeatable SPAR runstream without analytical gradients (RRSNG) 
is as foil ows : 


CXQT TAB 

TITLE"FUSELAGE MODEL 

TEXT 

••MEMBRANE ROD BEAM FUSELAGE MODEL 
••REPEATABLE PART WITHOUT GRADIENTS 

UPDATE*! 

$ MERGE NEW PROPERTIES HERE. 

UPOATE-0 
CXQT TQPO 
CXQT E 
CXQT EKS 
CXQT K 
CXQT INV 
CXQT M 

RESET G*QB1. 

CXQT AUS 

SYSVEC, ‘APPLIED FORCES 1 
1*1; J*6556Q;77; 10000.-20000. 20000. 
SYSVEC;UNIT VEC 
l»l; J=l»80; 1.0 
DEFINE WT-DEM DIAG 0 0 
DEFINE UN»UNIT VEC 
OBJFUN*XTY(UN,WT) 

CXQT OCU 
PRINT 1 OajFUN 
CXQT SSOL 
CXQT GSF 
CXQT PSF 
CXQT VPRT 

PRINT APPL FORC 1 1 
PRINT STAT DISP 1 1 


CXQT DCU 

PRINT 1 

STAT 

DISP 

1 

1 

PRINT 1 

STRS 

E21 

1 

1 

PRINT 1 

SIRS 

E23 

1 

1 

PRINT 1 

STRS 

E41 

1 

1 

TOC 1 

COPY 1#4 

OBJF 

AUS 

1 

1 

COPY 1,4 

STRS 

E21 

1 

1 

COPY 1,4 

STRS 

E23 

1 

1 

COPY 1,4 

STRS 

E41 

1 

1 

TOC 4 


CXQT EXIT 
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The End Processor 


The end processor, EPFUS, is also executed on the CDC 6600. The 
objective function and stress data saved on the SPAR library file, SPARED, are 
input to the end processor using FORTRAN-cal 1 able subroutines from SPAR (ref. 
13). Another file, ENDINXX, is stored on the CDC 6600 in NAMELIST form and 
input to the end processor. This file contains the number of each different 
type of element used in the model and the limits imposed on the behavior 
variables in constraint function evaluation. The format of the ENDINXX file 
is determined by the user because the user writes the end processor program 
dependent upon what elements are used in the finite-element model. The end 
processor then computes the constraints and outputs the objective function and 
constraints in a suitable format for input to the optimizer. The output is 
placed on file CNMNIO and transmitted to the PRIME on file TCMNIO. An example 
of an end processor program is listed in Appendix D. A sample of an ENDINXX 
file is as follows: 

$EPIN 

E23AL=2|300., 

E21AL=2000., 

E41AL=2000., 

NSE23=58, 

NSE21=76, 

NSE41=56, 

$END 

Gradients 

Options 1.3 and 2.3 have the capability for calculating gradients 
analytically. There are three problem independent programs used 
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to compute the analytical gradients. The first, BLDELDS, builds a file for 

input into the analysis program to calculate the derivatives for the mass and 

stiffness matrices with respect to the design variables in the nonrepeatable 

part. A listing of BLDELDS is in Appendix E. The program is stored on the 

CDC 6600 and is called by the TONOSNR (described in the section on the 

nonrepeatable analysis) sent to the CDC 6600 from the PRIME. The file 

containing the nonrepeatable SPAR runstream NRRS, is edited using the CDC Text 

Editor (ref. 12) commands from file EDITl to remove all but the element 

connectivity data. A listing of EDITl follows: 

F:/[XQT ELD/ 

F:/[XQT/;2 

D;* 

R 

F:/[XQT ELD/ 

E;* 

R 

D;* 

R 

ADD 

END 
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Another file, INPT, stored on the CDC 6600 and created by the user is 

input to BLDELDS along with the edited NRRS file. The INPT file contains two 

card types. The first type (only one card) has the format ((6(1X,I4) ,1X,A4) 

and contains the following variables: 

NOLC - number of load cases 
NODV - number of design variables 

ISNOLC - new starting number for load cases with derivatives 
JOINTS - number of joints in the model 
NDF - number of degrees of freedom 
NOEL - number of different elements 

VORB - a four character word to determine the type of analysis 
(e.g., STAT, VIBR, BUCK) 

The second type has the format (6(1X,A3,1X,I3,1X,I3)) consists of one or more 
cards, and contains the following variables (one set for each element): 

EL - element type (e.g., E21, E43) 

NSECT - largest section property number used for that element type 
NODVPE - number of design variables for that element type 

A sample listing of INPT follows. 

1 3 100 80 53 STAT 

E23 1 1 E21 1 1 E41 1 1 

BLDELDS creates a SPAR runstream and outputs it on file RSOUT. File RSOUT is 

then edited using the CDC Text Editor and commands from file EDIT2 to remove 

any extraneous blanks or characters before being input to SPAR. SPAR is then 

executed using RSOUT as input. The results are saved on a SPAR library file 

(NRLA) for later use in the repeatable part of the analysis. 



A listing of the EDIT2 file follows: 

RS:/ /,/ /;* 

RS: /,/ /;* 

A:/T0P0/;* 

>RESET MAXSUB=2500> 

END 

The other two gradient programs, GRDXXXX and DRVXXXX, are used in the 
repeatable part of the analysis. Listings of GRDXXXX and DRVXXXX are in 
Appendices F and G, respectively. Two procedure files SENDPRCG and TONOSG 
(similar to the SENDPRC and TONOS files previously described) are used to 
execute these programs in conjunction with the SPAR analysis program on the 
CDC 6600. The first part of the repeatable analysis with analytical gradients 
is identical to repeatable analysis described above. The CNMNIO file created 
by the optimizer is appended to the TONOSG file by SENDPRCG before it is 
transmitted to the CDC 6600. The GRDXXXX, DRVXXXX, INPTXX, CONSXX, and EDGRDS 
files are all stored on the CDC 6600. The contents of the INPTXX and CONSXX 
have already been described. Each of the GRDXXXX, DRVXXXX, INPTXX, and CONSXX 
general file names are replaced by specific names from the NAMEij file before 
execution. GRDXXXX uses the repeatable SPAR runstream (RRSG) as with 
analytical gradients to create a runstream to calculate the derivatives of the 
stiffness and the mass matrices with respect to the design variables when an 
element (such as a beam or plate) has more than one contributing factor. The 
program calls user supplied subroutine(s) with any of the following names 
DKDVE21, DKDVE22, DKDVE33, and/or DKDVE43. These subroutines compute the 
derivatives of the stiffness matrix with respect to a design variable for a 
particular element type in SPAR (ex. E21 or E43). The name(s) used depend 
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upon the elements used in the finite-element model. If a subroutine is not 
used, it remains as an unsatisfied external. Two integer parameters used in 
naming the created data sets are passed to each subroutine. The first is the 
counter for the design variable and the second is the number of unconstrained 
degrees of freedom (from 1 to 6). A listing of a sample DKDVE21 subroutine is 
shown in Appendix H. The subroutine card for DKDVE21 is as follows: 

SUBROUTINE DKDVE21(NDVJIM,NDF) 


where 


NDVJIM - the number of the design variables (first, second, third, etc.) 
NDF - number of degrees of freedom per joint squared 

A listing of the SENDPRCG file follows. 

SAVE TEMP 
DELETE UPDTE 
APPEND SENDFILE SPFPOUT 
APPEND SENDFILE CNMNIO 
DELETE CNMNIO 
RJSEND SENDFILE CDC 
DELETE SPFPOUT 
R TEMP 
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A listing of the TONOSG file follows. 

J08,T100j.CM10000Q » 

USER»USERNO> PASSWRD. 

CHAR6E>CHARGEND#LRC. 

' DELIVER. $UFO>SUBUFO>TCMNIO 
RFL, 100000. 

REDUCE(») 

MAP j OFF. 

COPYCR, INPUT, SPFPOUT. 

REWIND, SPFPOUT. 

GET,MERGFP/UN=753^37N. 

GBT,RUNSTREAM. 

EDIT,RUNSTRE AM,,MERGFP,E0DUT. 

REWIND, RUNSTREAM. 

GET,SPAR=SPAR1AI,DCU=DCU1AI/UN=750756N. 

get,sparla»sparname. 

SPAR,RUNSTRE AM,SPAROUT. 
RETURN,RUNSTREAM,SPFPOUT,MERGFp,EDnUT. 

REWIND, SPARLA, SPARED. 

GET,SPARLI8/UNs319925N. 

GET,GRDXXXX,DRVKXXX, INPTXX,CONSXX. 
GET,EDGRDS/UN*753A37N. 

CQPYCP, INPUT, CNMNIO. 

REWIND, CNMNIQ. 

GRDXXXX, INPTXX,CQNSXX,CNMNIQ, RSO'JT. 

REWIND, RSOUT. 

EDIT,RSnUT, , EDGRDS,EDQUT. 

REWIND, RSOUT . 

SPAR, RSOUT, SPAROUT. 

RETURN, EOGRDS, EDOUT, GRDXXXX . 
IFE,(NSUBS.NE.O),DERIV. 

REWIND,SPARLA,SPARLD,SPARLC,CNMNI0,C0NSXX, INPTXX. 
LDSEK LIBaSPARLI8,PRESET=ZER0) 

ORVXXXX, inptxx,consxx,cnmnio. 
RETURN,DRVXXXX,INPTXX,CONSXX, SPARLA,SPARLC. 
ENDIF,DERIV. 

REWIND,SPARLD, SPARLIB,CNMNID. 

GET,EP XXXX, ENDINXX. 

LDSET(LIB=SPARLIB,PRESET=ZERO ) 
EPXXXX,ENDINXX,CNMNIO,BLK. 

RETURN, SPARED, EPXX XX, ENDINXX, SPAR LIB. 

REPLACE, SPAROUT=OUMOUT,CNMNIOaDUMIO. V 

• REWIND, BLK,CNMNia. 

C0PYSBF,8LK, OUTPUT. 

COPYSBF,CNMNIO, OUTPUT. 

DAYFILE(L=DI$OUT) 

REPLACE, DISOUT. 

EXIT. 

DAYFILE(LaDlSOUT) 

REPLACE, DISOUT. 
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An example of a repeatable SPAR runstream with analytical gradients (RRSG) 
follows. 


CXQT TAB 

TITLE"FUSELAGE MODEL 
TEXT 

"MEMBRANE ROD BEAM FUSELAGE MODEL 
"REPEATABLE PART WITH GRADIENTS 
UPDATE-1 

S MERGE NEW PROPERTIES HERE. 

UPDATE=0 
CXQT E 
CXQT EKS 
CXQT TOPO 
CXQT K 
CXQT INV 
CXQT M 

RESET G=981. 

CXQT AUS 

SYSVEC; APPLIED FORCES 1 
1*1; J»65j69;77; lOOOO. -20000. 20000. 
SrSVEC;UNIT VEC 
I«lj 1.0 

DEFINE WT-DEM DIAG 00 
DEFINE UN*UNIT VEC 
OBJF AIJS 1 l = XTr(UN»wT) 

CXQT DCU 

PRINT 1 OBJF AUS 1 1 
CXQT SSOL 
CXQT GSF 
CXQT PSF 
CXQT VPRT 

PRINT APPL FORC 1 1 
PRINT STAT DISP 1 1 
CXQT OCU 

CHANGE l^STRS E21 1 1»FAMS E21 1 1 
TOC 1 • 

COPY 1#A OBJF AUS 11 
COPY 1,A FAMS E21 1 1 
COPY 1#A STRS E23 11 
COPY 1,A STRS EAl 1 1 



The SPAR runstream is output on file RSOUT. A file, EDGRDS, containing CDC 
Text Editor commands is used to remove any extraneous blanks or characters 


from RSOUT. SPAR is executed to compute the derivatives using RSOUT as 
input. The output data from SPAR, including stresses, stress derivatives, 
forces and moments, and derivatives of forces and moments are output on a SPAR 


library file, SPARED. A listing of EDGRDS follows. 


RS:/E+/,/+/;* 

RS:/. ./,/. ./;* 

RS:/. 0/,/. 0/;* 

RS:/. -/,/. -/;* 

RS:/ /./ /;* 

RS:/E-/,/-/;* 

RS:/ /,/ /;* 

RS:/ /./ /;* 

RS:/ W /,/ W/;* 

RS:/ G /,/ G/;* 

RS:/ W /J,W/;* 

RS:/ F /,/ F/;* 

RS:/ L /,/ L/;* 

END 


Program DRVXXXX is called if there is a need to convert forces and 
moments, and derivatives of forces and moments to stresses and stress 


derivatives. The program computes stresses and stress derivatives using a 
user supplied subroutine named BMSTRS (for beams) or PLTSTRS (for plates). As 
before, if a name is not used it remains an unsatified external. Four integer 
parameters are passed to BMSTRS. They are (1) a switch, (2) a counter so 
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certain computations can be skipped if they are not needed, (3) a block 
counter for accessing an array, and (4) another switch to determine if the 
beam is the contributing factor to the stress derivative. The first three 
parameters are also passed to PLTSTRS. A listing of a sample BMSTRS 
subroutine is provided in Appendix I. The subroutine for BMSTRS is as 
follows: 

SUBROUTINE BMSTRS ( ISW,KCNT, JCNT, IBEAM) 

where 

ISW,KCNT - a switch and a counter used to store certain beam data (ex. 
moments of inertia) and make certain computations (ex. 
derivatives of Y with respect to the design variables) only 
on the first time BMSTRS is called from DRVXXXX 

JCNT - a counter to find the location in blank common for various 
data times, depends on the number of beam elements 

IBEAM - 1 if beam is a contributing factor to stress derivatives 
- 0 otherwise 

This program also uses FORTRAN-cal lable SPAR subroutines (ref. 13) to input 
and output this data to and from the SPAR library, SPARED. After DRVXXXX 
completes execution the end processor is executed before control is returned 
to the PRIME. 

EXAMINING THE INTERMEDIATE RESULTS 

One of the key features in the distributed PROSSS is the capability for 
the user to examine the intermediate results. Because of this feature, the 
user does not have to wait until the end of the complete optimization process 
to determine whether or not the model, variables, and/or process are behaving 
as expected. The user, after examining the intermediate results, can let the 
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process continue; stop the process; or temporarily stop the process, make some 
changes to the optimization parameters or design variables, and restart the 
proccess. This gives the user much more control over the flow of the problem 
and should significantly reduce the total time required to reach the optimal 
solution. 

All examining of intermediate results is done on the PRIME because of the 
faster transmission rate (9600 BAUD) to the interactive terminals. Two files, 
CONOUTHOLD, and PLTDATA, created by the CONMIN driver program are essential 
for examining these results. The CONOUTHOLD contains a cumulative listing of 
all information output from CONMIN. This information is originally output on 
file CONOUT. At the end of the CNMNPRCl procedure file, CONOUT is appended to 
CONOUTHOLD to create the cumulative listing. It is best to examine the 
CONOUTHOLD file with the PRIME editor (ref. 14). The user is usually 
interested in seeing te progress of the objective function and design 
variables. Also of interest are the active and violated constraints. It is 
useful to check for errors and warning messages issued by CONMIN so that the 
process can be stopped at this point. 

If the user desires to see a graphical display of the optimization data, 
the PLTDATA file should be used in conjunction with program plots. The program 
plots either the objective function, design variables, or constraints versus 
the iteration number. The capability of listing these values in tabular form 
is also available. THe user makes a choice from a menu displayed on the 
screen (figure 3). The user can also decide whether or not the objective 
function and design variables are to be normalized. Examples of each type of 
plot are shown in figures 4-6. 
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RESTARTING PROSSS 


The user may wish to restart PROSSS for a particular point in the process 
because the computers go down, there is an error in the system, or the 
examination of the intermediate results revealed a change in the system, or 
the examination of the intermediate results revealed a change is needed in the 
model or optimization parameters. To restart PROSSS is quite simple. All of 
the necessary files required for restarting are saved by the procedure file 
CNMNPRCl. The user edits the first line of the NAMEij file to restart. The 
NRSTRT variable is set to 1. The NAMES variable is changed to create a unique 
name. The COUNT variable may or may not be updated depending upon whether or 
not the user has need of keeping track of the number of iterations through the 
system. Other names may be changed in the NAMEij file as needed. After these 
changes have been made, the user restarts PROSSS. After initialization, 
processing begins at CONMIN. 

SAMPLE PROBLEM 

Each option was executed to determine the final objective function 
(minimum mass, kg) of the finite-element model of the fuselage shown in figure 
7 using the NAMEij input files. The model is composed of 80 joints, 58 rods, 
76 beams, and 56 membranes. There are 352 degrees of freedom. The three 
design variables are (1) the cross-sectional area of the transverse stringers 
(beams), (2) the cross-sectional area of the longitudinal stringers (rods), 
and (3) the thickness of the panels (membranes). The design variables are 
handled by the optimizer in reciprocal form to improve the convergence. Their 
initial reciprocal starting values are 0.05 cm2, i,o cm2, 4.0 cm, 
respectively. The initial objective function is 5460.12 kg. Figure 8 shows 
the results for each option. All final objective functions are within 
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reasonable limits, even though some of the design variables differ 
significantly. The piecewise linear approach is 50 to 75 percent less 
expensive to execute than the nonlinear approach. Reference 6 describes, in 
detail, the results of this test compared to the results from an identical 
model tested on the mainframe-only version of PROSSS. The comparison is done 
in terms of accuracy, cost, CPU time, and total time. 
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GRADIENTS COMPUTED 


IN OPTIMIZER BY 

EXTERNAL TO OPTIMIZER BY 

METHOD 

FINITE DIFFERENCE 

FINITE DIFFERENCE 

ANALYTICAL 

NONLINEAR 

PROGRAMMING 

1.1 

1.2 

1.3 

PIECEWISE 

LINEAR 

PROGRAMMING 

N/A 

2.2 

2.3 


TABLE I.- OPTIONS FOR OPTIMIZATION. 
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FIGURE 2.- FLOl^CHART OF A TYPICAL OPTIMIZATION OPTION. 














ir^PUT NUMBER CORRESPONDING TO TYPE OF PLOT DESIRED 

1 - PLOT ITERATIONS VS OBJECTIVE FUNCTION 

2 - PLOT ITERATIONS VS DESIGN VARIABLES 

3 - PLOT ITERATIONS VS CONSTRANTS 
^ - PRINT NUMERICAL RESULTS 

5 - END 


FIGURE 3.- FENU FOR PLOTTING PROGRAM. 




FIGURE 4,- CHANGE IN OBJECTIVE FUNCTION FOR THE FUSELAGE MODEL. 
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FIGURE 5." CHANGE IN DESIGN VARIABLES (RECIPROCALS) FOR THE FUSELAGE MODEL. 
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FIGURE 6.- ACTIVE AND VIOLATED CONSTRAINTS FOR THE FUSELAGE MODEL. 





FIGURE 7.- 352 DOF FUSELAGE mDEL USED FOR TESTING. 
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wm 

D 
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D 

M/0 

D 

1.1 

1.6332 

1.6292 

4.7341 

5.0582 

.1000 
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1.2 

1.5713 

1.5681 

8.4175 

8.3682 

.1000 

.1000 

1.3 

1.5f|58 

1.5175 

10.6723 

11.0132 

.1042 

.1083 

2.2 

1.5659 

1.5149 

1.1770 

8.7796 

.1746 

.2158 

2.3 

1.6215 

1.6134 

1.6722 

2.8777 

.1000 

.1000 


M/0 = MAINFRAME-ONLY SYSTEM D = DISTRIBUTED SYSTEM 


FIGURE 8.- COMPARISON OF DESIGN VARIABLES. 
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APPENDIX A.- LISTING OF OPTION 1.1 DRIVER PROGRAM AND GENERAL SUBROUTINES 


C DISTRIBUTED PROSSS OPTll 
C 

INTEGER^^f COUNTiNAMES 
INTEGER+2 STAT#SW#ERROR> NRSTRT 
SW » 0 

SET UP INITIALIZATION PROCEDURE 

CALL SRCH$S(1,*NAME11'»6»1^0>IC00E ) 

READ! 5# 2) NRSTRT# ICR> NR »NOPT> NAMES* COUNT 
2 F0RMAT<3(I1,1X)* I2,1X*A4*1X*I3> 

IF (NRSTRT.EQ.O) CALL I NI T( NR* ICR * NAMES* NOPT ) 
COUNT » COUNT-1 

SET UP FOR RESTART IF NRSTRT » 1 

IF(NRSTRT.EQ.O» GO TO A 
CALL SRCH$$(A*«NAME11»*6*1*0*ICTOE) 

CALL RSTRT(NOPT*NAMES) 

CALL COMI$$( >PSTRTPRC'*8*6* ICQOE) 

CALL EXIT 

CALL SRCH$$(A, •RSTRTPRC»*8,1*0*IC0DE) 

CALL SRCH$$(5*«RSTRTPRC'*8*1*0*IC0DE) 

GO TO 6 

INITIALIZATION PROCEDURE FILE 


4 CALL SRCH$$( A, •NAM£11**6*1*0*ICQDE) 
CALL COMI$${ 'INITPRC**7*5*IC0DE) 

CALL EXIT 

CALL SRCH$$(A, • IN ITPRC * * 7* 1 #0* ICODE ) 
CALL SRCH$$(5* • INITPRC* *7*1*0*IC00E) 

SKIP IF NONREPEATABLE PART NOT REQUIRED 

IF(NR.EQ.O) GO TO 5 

CHECK STATUS OF RETURNED FILE 

CALL STCHEK(STAT*1) 

CALL SRCH$${5* »TNREPT**6*1*0*IC0DE J 
IF(STAT.EQ.O) GO TO 5 
ERROR « 2 
GO TO AO 
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GONMIN PROCEDURE FILE 

5 COUNT * COUNT+1 

IF( C0UNT.lt. 999) GO TO 6 
ERROR » 1 
60 TO 40 

6 CALL C0MI$$( 'CNMNPRCIS 8>6, ICODE ) 

CALL EXIT 

FRONT PROCESSOR PROCEDURE FILE 

CALL C0MI$S( »TFRTPRCS7»6# ICGDE) 

CALL EXIT 

SUBROUTINE TO UPDATE SEND FILE NAME (UNI3UELY) 

CALL UPDSF(CQUNT>SW»INFO) 

SW s SW+1 

UPDATE NAME PROCEDURE FILE 

CALL C0MI$$( «UPDTES5,6> ICODE) 

CALL EXIT 

SEND TO NOS PROCEDURE FILE 

CALL C0MI$$( 'TSENOPRCSB^b, ICODE) 

CALL EXIT 

CHECK STATUS OF RETURNED FILE 


CALL STCHEK(STATjO) 

IF(STAT.GT.O) GO TO 35 
ERROR » 2 
GO TO 40 

TEST FOR SYSTEM TERMINATION 

35 CALL SRCH$$(6# •60N0G0»,6,1,0^IC0DE ) 
IFUCODE.EQ.O) GO TO 100 

EDIT CNMNIO FILE RETURNED FROM NOS PROCEDURE FILE 

CALL C0MI$$( »EDPRC*>5#6^ ICODE) 

CALL EXIT 
GO TO 5 

PRINT ERROR MESSAGE AND STOP PROGRAM 

40 CALL ERRMSG(ERROR) 

100 STOP 
END 
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SUBROUTINE INIT ( NR » IC R > NAME 1 » NOPT ) 

CREATES A PROCEDURE FILE (INITPRC) FOR 
INTIALIZING files 

INTEGER* A NAME> NAME1> NAME2» STARTX> CONPAR# UFD> SUBUFD 
INTEGER*A S PL A# US ERNO# P ASS W # CHAR GE> INPT#C0NS 
DIMENSION NAME(2)#STARTX{2)>C0NPAR(2)>UFD(2),SUBUFD(2)>SPLA(2) 
DIMENSION USERN0(2)»CHAR6E(2)»INPT(2)>CQNS(2) 

OPEN PROCEDURE FILE 

CALL SRCH$$(2#*INITPRC*,7»2,0,IC0DE) 

WRITE(6>7) 

7 FORMATCSAVE TEMP'/*0PEN CONOUTHOLD 1 2»/»C 1*) 

SET UP TEMPORARY FILES FOR CNMNIO^ CONREST# AND BLK 

IF ( NOPT.EQ.il. OR. NOPT.EQ. 12. OR, NOPT.EQ. 13) WRITE(6> 8) 

0 FORMAT! »ED'/ ‘X* //*FIL C NMNI 0* / » E D * / » X • // • F IL CONREST* ) 

IF(NOPT. EQ. 22. OR. NOPT.EQ. 23 ) WRITE(6>9) 

9 FORMAT('ED»/*X*//*FIL BLK») 

INITIALIZE PASS» GONOGO# AND CNT 
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WRITE(6>10) 

10 F0RMAT(*ED*/ 'IV/'FIL PASS') 

IF<N0PT.EQ.12.0R.NQPT.EQ.a2 .0R.N0PT.EQ.23 ) WRITE (6^ 16) 

16 FORMAK '£0* / * l.OE 00 l.OE 00 l.OE 00 0«5E-01 * 

1 //'FIL CNT* ) 

INITIALIZE SENDPRC 

IF(N0PT*EQ. 11.0R.NQPT.EQ.12.QR.N0PT.EQ. 22) WRITE(6^ A5) 

^5 FORMATCED SENDPRC) 

IF(N0PT.EQ.13.0R.N0PT.EQ.23) WRITE(6> 46) 

46 FORMATCED SENDPRCG* ) 

WRITE(6/61) NAMEl 

61 F0RMAT( •C/SENDFILE/',A4, */100'/«FIL TSENDPRC) 

STORE NAMES FOR STARTX AND CQNPAR 



REA0(5,30) 

STARTX 


30 

FQRMAT(2A4) 

MRITE(6,40) 

STARTX 


40 

FORMATCED 

•>2A4/CFIL 

STARTX* ) 


READ(5>30) 

CONPAR 



WRITE(6,43) 

CONPAR 
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FORMATCED 

»>2A4/CFIL 

CONPAR * ) 


INITIALIZE front PROCESSOR 

READ(5»30) CONS 
WRITE<6#44) CONS 

44 FORMATCED FRONTPRC • / ' C /C QNSXX/ • , 2 A4C / 100 ' / • F I L TFRTPRC) 
INTIALIZE TONOS 
READ(5»30) NAME 

IF(N0PT.EQ.11.DR.N0PT.EQ.12.0R.N0PT.EQ.22) WRITE(6#65) NAME 
65 FORMATCED TONOS ' / ' C /RUNSTR E AM/ ' » 2 A4C /lOO * / * T • ) 

IF (N0PT.EQ.23) WRITE(6»64) NAME 
64 FORMATCED TONOSG • / • C /RUNST RE AM/ • * 2A4C /lOO • / « T* ) 
IF(N0PT.NE.13) GO TO 644 
WRITE(6>640) NAME 

640 FORMATCED TONOSl 3» / • C/ RSGR AD/ N 2A4, • / 100 ' /* T * ) 

READ(5,30) NAME 

WRITE(6»641) NAME 

641 FQRMATCC/RSNGRAD/'^2A4C/100'/*T» ) 

READ< 5>30) NAME ' 

WRITE(6»642) NAME 

64 2 F0RMATCC/EPYYYY/',2A4, */100«/'T* ) 

644 READ(5,30) UFD 

WRITE(6,66) UFD 
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66 FORMAT( 'C/UFD/'»2A4»*/100«/*T') 

REA0(5j30) SUBUFD 
WRITE<6,67) SUBUFD 

67 FORMAK 'C/SUBUFD/' »2AA> • /100*/»T* ) 

READ(5»30) NAME 

WRITE(6#68) NAME 

68 FORMAT ('C/DUMOUT/S2A A, '/lOO’WT* ) 

READ(5,30) NAME 

WRITE(6,69) NAME 

69 F0RMAT(*C/DUMia/'»2A4»»/100»#/«T* ) 

REA0(5»30) SPLA 

WRITE(6>70) SPLA 

70 F0RMAT(*C/SPARNAME/'»2A4>'/100»/*T* ) 

MRITE{6>71) NAMEl 

71 FORMAT { »C/J0B/*#AA,*/100'/*T* ) 

READ(5,30) USERNO 
WRITE(6^72) USERNO 

72 FQRMAT(*C/USERN0»/'»2A4» '/lOO'/'T* ) 

READ(5>30) PASSW 

WRITE(6#720) PASSW 

720 FORMAT( 'C/PASSWRD/S A4, »/100*/»T* ) 

READ(5>30) CHARGE 
WRITE(6»73) CHARGE 

73 F0RMAT(*C/>CHARGEN0»/'#2AA> »/100*/'T' ) 

READ(5>30) NAME 

WRITE(6,75) NAME 

75 FORMAT( ‘C/EPXXXX/S2AA# ' /lOO' W*T» ) 

READ(5#30) NAME 
WRITE(6»7A) NAME 

7A FORMAT( 'C/ENDINXX/'#2AA, '/lOOS/ »T» ) 

SKIP IF NO GRADIENTS ARE REQUIRED 

IF (NOPT.EQ.n.OR.NOPT.EQ.12*QR,NOPT.EQ.22) 60 TO 755 
REA0(5»30) NAME 
WRITE(6#750) NAME 

750 FORMAT! ‘C/DRVXXXX/ *#2A4, */100'/'T» ) 

READ(5»30> NAME 

WRITE(6>751) NAME 

751 FORMAT! *C/GRDXXXX/S2AA, • /100‘/'T* ) 

READ!5,79) NSUBS 

79 F0RMAT!I1) 

WRITE(6#752) NSUBS 

752 FORMAT !‘C /NSUBS/'* I 1»*/100»/‘T») 

READ!5»30) INPT 
WRITE!6*753) INPT 

753 FORMAT! • C / INPTXX/ • * 2 A A, • /lOO • / *T • ) 

WRITE!6*754) CONS 

75A F0RMAT!'C/C0NSXX/'»2A4, »/100'/'T' » 
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REMOVE EXTRANEOUS BLANKS 
755 WRITE(6,76) 

76 F0RMAT<»C/ / / 100 • / ' T ' / ' C / //lOO* /' T' /* C/ //lOO*) 
MRITE(6#77) NAMEl 

77 FORMAK *T •/ 'C/ / / 100 ' / • T ' / « C / //100*/*FIL * # AA ) 

INITIALIZE STATIN 

WRITE(6#78) NAMEl 

78 F0RMAT(«ED'/A4//'FIL STATIN*) 

IF NOCR EQ 0 JUST CHANGE NAMES FOR PROGRAMS 
IF NOCR NE 0 EDIT FILES TO COMPILE PROGRAMS 


FRONT PROCESSOR 

READ(5#115) NAME#N0CR 
115 F0RMAT(2AA,1X^ ID 

IF(NOCR.NE.O) WRITE(6,80) NAME 

80 FORMATTED CRF PROC • / • C /F PXX XX/ • , 2 A AD / 1234 • / • F I L TCRFPROC*) 
IF(NOCR.NE.O) WRITE(6»120) 

120 FORMAT(»CO TCRFPROC 6*/»C 6*/*DELETE TCRFPROC*) 
IF(NOCR.EQ.O) WRITE(6,800) NAME 
800 FORMAT(*CN ♦»#2AAD ♦F PR OC E S SOR • ) 

EVAL 

IF(N0PT.EQ.11.0R.N0PT.EQ.13) GO TO 860 
READ(5#115) NAME#N0CR 
IF(NOCR.NE.O) WRITE(6D5) NAME 

85 FORMATCED C RE VAL • / ' C / EVAL XX/ ' >2 A A# ' / 1234 * / * F IL TCREVAL*) 
IF( (NOPT. EQ.12.0R.N0PT.EQ.22.0R.N0PT.EQ.23) .AND.NOCR.NE.O ) 

1 WRITE(6,1A0) 

140 FORMATCC 6*/'C0 TCREVAL*/'C 6*/*DELETE TCREVAL*) 
IF(NOCR.EQ.O) WRITE(6>850) NAME 
850 FORMAT(*CN *#2A4D EVAL*) 

CONMIN 

860 READ15#115) NAME#NOCR 

IF(NOCR.NE.O) WRITE(6>90) NAME 
IF (NQCR.NE.O) WRITE(6,130) 

130 FORMAT(*C 6*/*C0 TCRCONMIN */*C 6*/*DELETE TCRCONMIN*) 
IF(NOCR.EQ.O) WRITE<6#900) NAME 
900 FORMATCCN *»2A4>* CONMIN*) 
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IF(N0PT.NE.12) GQ TO 960 


C RERITES 

e 

READ(5#115) NAME>N0CR 
IF(NOCR.NE.O) WRITE(6>950) 

?50 FORMATCCO CRRITE 6»/»C 6* ) 

GO TO 1A5 

SELECTS 

READ(5#115) NAME^NOCR 
IF(NOCR.NE.O) WRITE(6,951) 

951 F0RMAT(*C0 CRSEL 6*/»C 6M 

TEST TO DETERMINE IF MODEL CREATED WITH DATA GENERATOR 

1A5 IF(ICR.NE.O) CALL CRSNPT 

TEST TO DETERMINE IF NONREPEAT ABLE PART IS TO BE DONE 

IF(NR.NE.O) CALL NONRPT ( NR, S PLA» INPT> USERNQ»P ASSW» CHARGE* 
1 UFD*SUBUFD,NDPT) 

WRITE(6,150) 

150 FGRMATCR TEMP' ) 

CALL SRCH$$(A,«INITPRC',7,2,0, ICODE) 

RETURN 
END 

SUBROUTINE RSTRT(NOPT, NAMES ) 

THIS SUBROUTINE SETS UP A PROCEDURE FILE FOR 
RESTARTING THE PROGRAM 

INTEGER^A NAMES, NAMERS* NAMOLD 
DIMENSION NAMERS(2) 

CALL SRCH$$(2, 'RSTRTPRC • »8,2,0,IC0DE) 

WRITE(6,10) 

10 FORMATCSAWE TEMP*) 

IF(N0PT.EQ.11.0R.N0PT.EQ.12.0R.N0PT.EQ.13) GO TO 20 
IF(N0PT.EQ.22.0R.N0PT.EQ.23 ) 60 TO 60 

SET UP FILES FOR OPTIONS 1.1* 1.2* AND 1.3 

20 WRITE(6,30) 

30 FORMAT! 'DELETE PASS'/'CN SPASS PASS*) 

WRITE(6*A0) 

AO FORMATCDELETE CNMNI0'/*CN SCNMNIO CNMNIO* ) 

WRITE(6*50) 

50 FORMAT! 'DELETE CONREST'/*CN SCONREST CONREST*) 
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GO TO 100 


SET UP FILES FOR OPTIONS 2.2 AND 2.3 
60 WRITE(6#70) 

70 FORMAT( ‘DELETE CNT’/'CN SCNT CNT») 

WRITE(6>80) 

80 FORMAT(»DELETE BLK*/»CN SBLK BLK*) 

WRITE(6»90) 

90 FORMATC ‘DELETE STARTX‘/‘CN SSTARTX STARTX') 

REINITIALIZE SENDPRC 

100 WRITE(6»101) 

101 FDRMAT(»DELETE TSENDPRC) 
IF(N0PT.EQ.11.0R.N0PT.EQ.12.0R.N0PT.EQ.22) WRITE<6> 102) 

102 FORMAT(‘ED SENDPRC‘) 

IF (N0PT.EQ.13.0R .N0PT.E0.23 ) WRITE (6# 103) 

103 FORMAT(‘ED SENDPRCG‘) 

WRITE(6>107) NAMES 

107 F0RMAT<*C/SENDFILE/*»A4> »/100*/‘FIL TSENDPRC* ) 

REINITIALIZE STATIN 

CALL SRCH$$ ( 1, ‘ STATIN ‘,6#1>0,ICQDE) 

READ(5#110) NAMOLD 
110 FDRMAT(AA) 

CALL SRCH$$(4#‘STATIN‘#6^1>0>IC0DE ) 

MRITE(6»115) NAMOLDfNAMES 
115 FORMATCCN *,A4#1X»A4) 

WRITE(6»120) NAMES 

120 FORMAT(‘DELETE STATIN ‘ /' ED • /A4// ‘ F IL STATIN*) 

WRITE(6#130) 

130 FORMAT(‘R TEMP* ) 

CALL SRCH$$(4# ‘RSTRTPRC* ,8, 2,0, ICODE) 

RETURN 

END 

SUBROUTINE NONR PT ( NR , SP L A, I NPT,US ERNO, P ASS W,CHARGE, UFD,SUBUFD, 

1 NOPT) 

CREATES PROCEDURE FILE TO SEND TO CYBER FOR NONR EP EAT ABLE PART 
OF ANALYSIS 

INTEGER*4 NAME, INPT, SPL A,NR OPT, NRR S, USERNO, CHARGE, UFO, SUBUFO, PASSW 
1 ,TNREPT 

DIMENSION NRRS(2),SPLA(2), INPT(2),NAME(2),USERN0<2),CHARGE(2) 

1 ,UFD(2),SUBUFD(2) 

CALL SRCH$$(1, ‘NAME NR *,6,3,0, ICODE) 

WRITE<6,10) 
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10 F0RMAT(*ED TDNOSNRM 
READ(5>20) NRRS 

20 F0RMAT(2A4) 

WRITE(6#30) NRRS 

30 F0RMAT(*C/NRRS/*»2AA> */100'/*T') 

WRITE(6»50) NDPT 

50 F0RMAT(*C/NR0PT/*#I2#»/100* /»TM 

WRITE(6,60) INPT 

60 FORNAK 'C/INPTXX/'»2A<f, • /100‘ / »T* ) 

REA0(5#20) NAME 
WRITE<6#70) NAME 

70 FORMAK *C/FLX/S2A^# »/100*/ *T* ) 

HRITE(6>71) USERNO 

71 F0RMAT(*C/USERN0,/'»2AA, VIOO*/'!* ) 

WRITE(6>72) PASSW 

72 F0RMAT(*C/PASSWRD/‘>AA>»/100'/'T» ) 

WRITE(6,73) CHARGE 

73 FORM AT ( 'C /»CHARGENQ, / • , 2AA> '/lOO'/'TM 
WRITEC6#7A) UFO 

7A FORMAT(»C/UF0/'»2AA,»/lOO»/'T') 

WRITEC6#75) SUBUFD 

75 format ( 'C/SUBUFD/'» 2AA» '/lOO'/'T* ) 

READ(5>76) TNREPT 

76 FORMAT(AA) 

WRITE(6>77) TNREPT 

77 FORMA T( 'C /JOB/ ' »A4» • / 100*/ * T » ) 

WRITE(6^80) SPLA 

80 FORMAT ( • C /NRLA/ « , 2AA, '/lOO* /*T* ) 

WRITE(6#8A) 

8A F0RMAT(*C/ //100G'/*T*/*C/ //lOOG • / *T • / *C / //lOOG*) 
WRITE(6»85) TNREPT 

85 FORMATCFIL *#AA) 

WRITE<6»86) TNREPT 

86 FORMAT( 'ED*/AA//'FIL NRSTAT*) 

WRITE(6»90) TNREPT 

90 FORMAT( ‘APPEND '^AA#' TS PAR IN •/* DELETE TSPARIN*) 

WRITE(6,110) TNREPT 
110 FORMAT( *R JSEND *#AA#' CDC) 

CALL SRCH$$(A#*NAMENR*»6>3»0>IC0DE) 

RETURN 
END 

SUBROUTINE CRSNPT 

EDITS FILE CREATED BY DATA GENERATOR AND LOADS IN FILE CONTAINING 
PROBLEM DEPENDENT DATA 

DIMENSION FNAME<2),LNAME(2) 

INTEGER+A NDV>LNAME#FNAME 
CALL SRCH$$( 1»*NAMECRRS »»8f 3>0»IC0DE ) 
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READ IN NUMBER OF DESIGN VARIABLES 

READ(5#1) NDV 

1 F0RMAT(I3) 

READ IN NAME OF FILE CREATED BY DATA GENERATOR 

READ15»2) FNAME 

2 F0RMAT(2A4) 

WRITE(6,3) FNAME 

3 FORMATCED '»2AA/»'L JL DC ' / 'N ' / ' D • / ' L CMDM 
WRITE(6>4) 

4 FORMATC *D«/ *L CMD • / • C /C MD/ C XQT/ ' / » C>5 /13 > ' // • $ST ART* / / ) 

NDVMl = NDV-1 

IF(NDVMl.EQ.O) GO TO 8 
DO 7 I « IfNDVMl 
WRITE(6>6) 

6 FORMAT(»L E • / • N-1 • // ' SEND' / • SSTART* // ) 

7 CONTINUE 

8 WRITE(6#9) 

9 FORMATCL CMD ' / ' D • / • N-1 ' / / • $ END • / • C XQT EXIT'//'T') 

READ IN FILE NAME CONTAINING PROBLEM DEPENDENT DATA TO BE LOADED 
IN AT THE TOP OF THE FILE 

READ(5#2) LNAME 
WRITE(6/10) LNAME 

10 FORMATCLOAD ',2A4/,'FIL TSPARIN*) 

CALL SRCH$i(4, 'NAMECRRS* »8>3»0>IC0DE) 

RETURN 

END 

SUBROUTINE UPDS F ( CNT> SW> IN FO ) 

THIS SUBROUTINE CREATES A PROCEDURE FILE (UPDTE) 

TO CREATE UNIQUE JOB NAMES TO SEND TO NOS 

INTEGER+2 NME(2)#SW 
INTEGER*^ CNT,CNTM1 

MODIFY STATIN 

CALL SRCH$$(1»» STATIN',6#1»0»IC0DE ) 

READ(5^10) NME 
10 F0RMAT(2A2) 

CALL SRCHS$(4, ' STATIN* > 6, 1, 0, ICODE) 

CALL SRCH$S(2> 'UPDTE'»5,2»0»IC0DE ) 

WRITE(6,20) NME^NMEKNT 

20 FORMATCSAVE TEMP'/'ED STAT IN • / • N • / • C/ • > 2A2» • / « , 2A2» 13 ) 
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WRITE(6>30) 

30 FORMATCC/ //»/»C/ //»/'FIL TSTATIN') 

CNTMl ■ CNT-1 

CHECK FOR FIRST TIME THROUGH 
IF(SW.EQ.O) 60 TO 41 

MODIFY SENOPRC DEPENDING UPON PASS COUNTER 

IF(CNT.LT.IO) WRITE(6>35) NME > CNTMl^NME# C NT 

35 F0RMAT(»ED TSENDPRC • / • C / S 2 A2, II , • / * , 2A2» 1 1, • /lOO • / » T* ) 
IF(CNT.E0,10) WRITE(6#360) NME#CNTM1,NME> CNT 

360 FORMAT<»ED TSENDPRC • / » C / • » 2 A2/ Il> • / S 2A2> 1 2/ ' /lOO ' / *T • ) 

IF (CNT.GT. 10. AND.CNT.lt. 100 ) WRITE (6, 36) NME»CNTM1# NME>CNT 

36 FORMATCED TSENDPRC * / ' C / ' » 2 A2, I2> • / • » 2A2# 1 2> • /lOO ' / • T * ) 

IF(CNT. 6E.100.AND.CNT.lt. 1000) WRITE(6#37) NME# CNTM1» NME > CNT 

37 FORMAT(»ED TSENDPRC • / * C / * , 2 A2. 13, • / • , 2A2> 1 3# • /lOO • / * T • ) 
WRITE(6,38) NME, CNTMl 

38 FORMATCL DE L •//• DELETE »>2A2,I3//) 

IF(SW.NE.l) WRITE(6,39) 

39 FORMAT! 'NWD'/'T* ) 

IF(SW.EQ.l) WRITE(6,40) 

40 FORMATCT*) 

GO TO 49 

41 WRITE(6,45)NME,NME,CNT 

45 FORMAT!'ED TSENDPRC • / • C / ' , 2 A2, • / • , 2A2» 1 3, • /lOO* / » T • ) 

49 WRITE(6,50) NME,NME,NME,NME 

50 FORMAT! »C/*2A2, » / • 2 A2 , • / 99 ' /* T • / • C/ » 2 A2, • / • 2 A2, • 799 » , / • T* ) 

WRITE!6>60) NME,NME 

60 FORMAT! »C/» ,2A2, • / ' , 2A2, • 799* /• F IL TSENDPRC*) 

MODIFY FILE TO BE SENT TO NOS 

IF!CNT.LT.10) WRITE!6,70) NME , NME , NME, CNT, INFO, NME, CNT 
70 FORMAT!*ED ' 2A2, 7 • N • , 7 » C 7 • 2 A2, '7 • 2A2, 1 1, • 7 100 • 7 * T * 7 
1 •C7INF07',I1, '7100'7’FIL ',2A2,I1) 

IF (CNT. GE.10.AND.CNT.lt. 100 ) WRITE!6,80) NME,NME,NME, CNT, ' 

1 INFO,NME,CNT 

80 FORMAT(*ED ' 2A2, 7 ' N • , 7 » C 7 • 2 A2 , • 7 • 2 A2, 1 2, ‘ 7 100 • 7 ' T* 7 
1 •C7INF07',I1,»7100*7»FIL »,2A2,I2) 

IF!CNT. 6E.100.AND.CNT.lt. 1000) WRIT£!6,90) NME, NME, NME ,CNT, 

1 INFO, NME, CNT 

90 FORMAT!*ED » 2A2, 7 »N • , 7 • C 7 • 2 A2, • 7 • 2 A2, 1 3, • 7 100 • 7 • T • 7 
1 •C7INF07»,I1, »7100'7«FIL *,2A2,I3) 

WRITE!6,100) 

100 FORMAT!*R TEMP* ) 

CALL SRCH$$! 4, 'UPDTE* ,5,2,0,IC0DE ) 

RETURN 
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END 

SUBROUTINE ERRMSG( ERROR ) 

CALL SRCH$$(2, 'ERRFILE • , 7, 2#0, IC ODE ) 

IF (ERROR. NE.l) GO TO 15 
WRITE(6»10) 

10 F0RMAT(» MORE THAN 999 ITERATIONS. STOP PROGRAM.* J 
GO TO 100 

15 IF(ERR0R.NE.2) GO TO 100 
WRITE (6»20J 

20 F0RMAT(» CANNOT FIND COMET. STOP PROGRAM.') 

GO TO 100 

100 CALL SRCH$$(4, 'ERRFILE*»7,2#0> ICODE) 

RETURN 

END 

SUBROUTINE STCHEK(STAT# ICHK ) 

THIS SUBROUTINE CHECKS STATUS UNTIL FILE IS RETURNED 
FROM NOS 

INTEGER *2 SN AME ( 3 ) ^ RNAME ( 3 ) > UFD ( 3 ) > STAT 
INTEGER+4 TIMERjCNTR 
TIMER » 120000 
CNTR = 0 

PUT JOB TO SLEEP FOR TWO MINUTES 
CALL SLEEP$(TIMER) 

DETERMINE FILE NAME FOR NONRE P E AT A8 L E PART 

IF(ICHK.EQ.O) GO TO 5 

CALL SRCH$$(1, *NRSTAT'»6»1,0> ICQDE ) 

READ(5>10) SNAME 

CALL SRCH$$(4» »NRSTAT'»6#1»0» ICODE) 

CALL SRCH$$(5» ' NRST AT •> 6> 1 . 0» ICODE ) 

GO TO 20 

DETERMINE FILE NAME FOR REPEATABLE PART 

5 CALL SRCH$$(1» 'TSTATIN'^TU^O# ICODE) 

REA0(5#10) SNAME 
10 F0RMAT(3A2) 

CALL SRCH$$(4, 'TSTATIN' >7,1,0# ICODE) 

CHECK STATUS UNTIL STAT NE 0 

20 CALL COMSTA(SNAME,RNAME,UFD,STAT> 

IF(STAT.NE.O) GO TO 30 
CALL SLEEPS (TIMER) 
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INCREASE SLEEP PERIOD EVERY 5 ITERATIONS 

CNTR ■ CNTR+1 
IF(CNTR.LT.5) GO TO 20 
CNTR » 0 
TIMER » TIMER*2 
GO TO 20 
30 RETURN 
END 

SUBROUTINE COMSTA ( SNAME>RNAME#UFD> STAT ) 

COMSTAT -COMET STATUS. RETURNS STATUS OF FILE TRANSMITTED 
TO THE HOST CDC COMPUTER. 

CALLING SEQQUENCE s 

SNAME -THREE WORD ARRAY IN CALLING PROGRAM CONTAING 

NAME OF FILE SENT TO HOST. LEFT JUSTIF lED, BLANK FILL. 

RNAME -THREE WORD ARRAY IN CALLING PROGRAM WHERE NAME 
OF FILE RETURNED FROM HOST IS STORED. LEFT 
JUSTIFIED»BLANK FILL. 

UFD -THREE WORD ARRAY IN CALLING PROGRAM WHERE UFD OF RETURNED FILE IS STORED 
THE CALLER'S LOGIN UF D> OT HERWI S E » COMET. 

6CHARACTERS LEFT JUSTFI ED> BLANK FILL. 

STAT -STATUS OF FILE SEND TO HOST. 

ONE WORD INTEGER IN CALLING PROjpRAM. 

-1 » ERROR HAS OCCURED»STATUS CANNOT BE DETERMINED. 

THE CONTENTS OF RNAME AND UFD ARE NOT VALID. 

0 - FILE HAS NOT BEEN RETURNED. 

THE CONTENTS OF RNAME AND UFD ARE NOT VALID. 

1 ■ FILE HAS BEEN RETURNED. 

THE CONTENTS OF RNAME AND UFD ARE VALID. 

INTEGER*2 SN AME ( 3 ) # RN AME ( 3 ) # UF D( 3 ) » R J DATA ( 30 ) # ERR 
INTEGER»2 STAT, CODE# RNW» COMUFD (3 )# BLANKS ( 3 ) 

INTEGER*^ POS,EOFPOS 
SINSERT SYSCOM>KEYS.F 

LOGICAL NAMEQV, FOUND 
DATA FOUND/. FALSE./ 

DATA COMUFD/'COMET'/ 

DATA BLANKS/' '/ 

C ATTACH TO UFD COMET 

CALL ATTACH(COMUFD,INTS( :100000> , BLANKS, 0, $900) 

CALL TNOUCATT TO COM', 10) 

C WRITE(1,777)STAT 
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777 FORMATMSTAT » SIS) 

C TRY TO OPEN STATUS FILE ON UNIT 12 
5 CALL 0PRDWT( »STATUSS12>CQDE) 

. IF(CODE*NE.O) GO TO 950 

C OPRDWT CYCLES UNTIL 'STATUS' IS UPDATED DURING RE-WRITE 

CALL TNOUCOPEN STATUS', 11) 

SET PQS TO EOF AND THEN 
REWIND SO WE CAN HAVE 
A VALID EOF POSITION 

CALL PRWF $$( KSPOSN, 12# L QC(0), 0,1 0000000, RNW, CODE) 

CALL PRHF$$(K$RPQS,12,L0C(0),0,P0S,RNW,C0DE) 

EOFPQS»POS 


CALL PRW F$$(K$P0SN,1 2, LQC(0 ), 0,-10000000, RNW, CODE) 
CALL PRWF$$(K$RP0S,12,L0C(0),0,P0S,RNW,C0DE) 


READ A STATUS ENTRY AND LOOK FOR MATCH WITH SNAME 
CALL PRWF$$(K$RP0S,12, LOC ( 0 ), 0, POS , RNW, CODE ) 
IF(CODE,NE.O) GO TO 100 
IF(POS.GE,EOFPOS) GO TO 100 
CALL I$AD07(12,RJDATA,30) 

IF(.NOT.(NAMEQV(SNAME(l),RJDATA(i)))) GO TO 10 
NAME FOUND, CHECK STATUS 
CALL TNOUCFOUND NAME*, 10) 

FOUNDS. TRUE. 

WRITE(1,44^) RJDATAOO) 

IF(RJDATA(30) .LE.0)G0 TO 800 
A FORMAT( 'RJDATA 30 = ',15) 

CHECK FILE NAME FOR T$ PREFIX, IF PREFIX NOT 
PRESENT, AN ERROR HAS OCCIJRED. 

IF(RJDATA(17) .NE. 'TS» )G0 TO 800 

CHECK FOR REPEATABLE NAME TCMNIO 

IF (RJDATA (17 ) .EQ.'TC . AND. R JOATA ( 1 8 ) . EQ . ' MN ' . AND . 
1 RJDATA( 19) .EQ. '10' ) GO TO 19 

CHECK FOR NONREPEATABLE NAME TNREPT 

IF(RJDATA(17) .EQ.'TN'. AND.R JDATA( 18) .EQ.'RE'.AND. 
1 RJDATA(19) .EQ. *PT' ) GO TO 19 
GO TO 800 

C FILE HAS BEEN RETURNED, COPY NAME TO RNAME. 

19 DO 20 1*1,3 

RNAME( I)»RJDATA( I+16) 
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20 CONTINUE 

C GET UFD STATUS 

IF(RJDATA(30).EQ.1>G0 TO 50 
DO 25 1-1,3 

UFD(I)-COMUFD(I) 

25 CONTINUE 

GO TO 60 

50 DO 55 1-1,3 

UFDd )«RJOATA( 1 + 6) 

55 CONTINUE 

60 CONTINUE 

C SET RETURN STATUS TO FOUND 
STAT * 1 
GO TO 10 
100 CONTINUE 

IF<FOUND) STAT-1 
CALL ATTACH(0,0,0,0,0) 

CALL SEARCH(A,0,12) 

RETURN 

800 CONTINUE 

C SNAME NOT FOUND, SET STATUS AND RETURN 
STAT » 0 
DO 810 1-1,3 

RNAMEiD-BLANKSd) 

UFD( I)-BLANKS( I) 

810 CONTINUE 

< CALL TNOU(«NAME NOT FOUND', lA) 

FOUND-, FALSE. 

GO TO 100 
900 CONTINUE 

C CANNOT FIND COMET, SET ERROR STATUS AND RETURN 

STAT « -1 
GO TO 100 
950 CONTINUE 

C CANNOT OPEN STATUS, IF BUSY,TRY AGAIN, 

CALL 6ETERR(ERR,1 ) 

IF(ERR.EQ. 'SI» )G0 TO 5 
C FILE NOT FOUND MEANS NOT RETURNED YET. 

CALL TNOU(»CAN"T OPEN STATUS', 17) 

GO TO 900 
END 

C OPRDWT->OPEN READ AND WAIT 

SUBROUTINE 0PRDWT(NAME,UNIT,RC0DE) 

C 

INTEGER+2 NAME(3),UNIT,TYPE,CQDE, I,RCODE 
C 

$INSERT SYSCOM>KEYS.F 
DATA 1/0/ 

RCODE-0 
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10 CALL SRCH$S(K$EXST>*r$0001*>6/l#0#C0DE) 

C 

C IF T$0001 DOES EXIST THEN STATUS IS BEING RE-WRITTEN 

C SO LOOP UNTIL THRU WITH UPDATE FUNCTION 

C 

IF(CODE.EQ.O) GO TO 999 
15 DO 200 I»l^25 /O'KEEP TRYING*/ 

CALL SRCH$$(K$READ>NAME#6#UNIT^TYPE#C0DE) 

IF(CQDE.EQ.O) RETURN /*SUCC ESSFULL*/ 
IF(C0DE.E0.4.0R.C00E.EQ.5) GO TO 100 
RC0DE*1 
RETURN 

100 CALL SLEEP$< INTL(5000) ) 

C 

C SLEEP FOR 5 SECONDS 

C 

200 CONTINUE 
RCODE-1 
RETURN 

999 CALL SLEEP$(INTL(5000)) 

GO TO 10 
END 
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APPENDIX B.- LISTING OF MAIN PROGRAM FOR CONMIN. 


PROGRAM C0NMS1(INPUT»0UTPUT,TAPE5»INPUT>TaPE6»0UTPUT> 
1TAPE7» TAPES, TAP E9»TAPE10» TAPE 11) 

NEW CONMIN SUBROUTINE 

COMMON/CNMNl/D£LFUN,DABFUN,FDCH, FOCHM,CT,CTMIN,CTL,CTLMIN, 
1ALPHAX,AB0BJ1,THETA,0BJ,N0V,NC0N,NSIDE,IPRINT,NFDG,NSCAL, 
2LINOB J,ITM AX, ITRM, ICNDI R, IGOTO,N AC, INFO, INFOG, ITER 
COMMON/CNMN2/RDUM( 50),IDUM( 25) 

COMMON X(20 ), VLB(20),VUB (20),G(200),SCAL(20),DF(20), 
1A(20,100),S(20),G1(200),G2(200),B(100, 100),C(100),ISC( 200), 
2IC(100),MS1{200) 

COMMON/CONS AV/RSAV (50), I SAV (25) 

DIMENSION KSC(200),SAVX(20) ,SAVG(200) 

INTEGER*2 ICODE 
C REAOdO, PASSAGE) 

READ(11,10) NPASS 
10 FORMAT(Il) 

C first PASS, NPASS»1, SUBSEQUENTLY NPASS-2. 

GO TQ( 100, 200), NPASS 
100 CONTINUE 
C READ(5,C0NPAR) 

RE AD (5,*) I PRINT, NOV, ITM AX, NCON, NFDG, NS I DE, ICND IR, 
1NSCAL,LIN0BJ, ITRM, FDCH,FDCHM, CT,CTMIN,CTL,CTL MIN, 
2THETA,PHI,DELFUN,DABFUN, JSC,N1,N2,N3,N4,N5, 
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3ALPHAX»ABG0Jl>IGQTa, (VLB(n»I»l»NDV)> (VUB(I)> I»1>N0V) 

JITER « 0 
MSC » 0 

IF(JSC.LT.O) MSC » 1 
DO 15 I - 1>NC0N 
ISC(I) « MSC 
15 CONTINUE 

IF ( JSC .EQ.O.OR. JSC.EQ.-999) GO TO 13 
MSC = 1 

IF(JSC.LT.O) MSC « 0 
IF(JSC.LT.O) JSC » JSC*“1 

CALL SRCH$$( INTS(1)» •ISCFILESINTS(7)>INTS(A) »INTS(0)> ICODE) 
REA0(8,12) (KSC (I ) »I»1# JSC ) 

12 FORMAT(IOI^) 

CALL SRCH$$(INTS(A),'ISCFILESINTS(7),INTS(4),INTS(0),IC0DE) 
DO 11 I » 1# JSC 
ITEMP » KSCd) 

ISC(ITEMP) = MSC 
11 CONTINUE 

13 REWIND 13 

WRITE(13/1A) JITER»NDV,NCON> ( I SC ( U , I» 1» NC ON ) 
lA F0RMAT(20I5) 

C READ(8,STARTX) 

CALL SRCH$$(INTS(1)> 'STARTX V INTS ( 6 ) , INTS ( A ) , INTS ( 0 ) > I CODE ) 
READ(8#30) NDV> (X( I)» I=1»NDV) 

30 F0RMAT(I5/(6E13.5) ) 

CALL SRCH$$(INTS(A), »STARTXSINTS(6)#INTS(4), INTS(0)»IC0DE) 

GO TO 201 
200 CONTINUE 
C READ(7#SAVE) 

RE AO (7 ) IPRINT>NDVdTMAX>NC0N,N3IDE»ICN0IR,NSCAL#NF0G> 

1FDCH» FOCHM# CT#CTMIN>CTL,CTLMIN>THETA>PHI>NAC#DELFUN>DABFUN» 
2LIN0BJ,ITRM» ITER» INFOG# I GOTO# INF 0, OB J # 

3RDUM»IDUM# 

AX#OF,G#ISC> IC» A,S>G1#G2»C>MS1»B# VLB # VUB# SC AL# RSA V# I SAV #NCOUNT 
5#N1#N2#N3#N4#N5#ALPHAX# ABQ8 J1 
REWIND 7 

C REA0(9#LINKE ) 

REA0(9#A0) OBJ#G 
AO F0RMAT(E13.5/(6E13.5) > 

ITERSV « ITER 
DO A1 I « 1»NC0N 
SAVGd) * G(I) 

A1 CONTINUE 

DO A2 I « 1#NDV 

SAvxn ) - x(i ) 

A2 CONTINUE 

SAVOBJ - OBJ 
SAVCT • CT 
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SAVCTL - CTL 
REWIND 13 

READ(13#14) JITER#NDV>NC0N» {ISC(I)>I»1#NGQN) 

JITER - JITER+1 
201 CONTINUE 
NPASS-2 
REWIND 11 

C WRITE(10#PASSAGE) 

WRITE(11>10) NPASS 
C SOLVE OPTIMIZATION 

CALL CONMIN(X> VLB» VU3»G»SCAL# DF» A^S»G1#G2# BfC# 
♦ISC#IC>MS1#N1,N2#N3#NA#N5) 

IF (JITER. EQ.O) GO TO A8 
IF(JITER.EQ.l) GO TO 47 
IF( ITER.EQ.ITERSV) GO TO 48 
47 ITERSV ■ ITER 
REWIND 13 

WRITE (13# 14) JITER#NOV#NCQN#( ISC (I )#I»1>NCDN) 

REWIND 14 

WRITE (14, 45) SAVOBJ, (SAVX( I ), I»l, NOV), SAVCT, SAVCTL 
1 , (SAVG(I),I»1,NC0N) 

45 F0RMAT(6E13.5) 

C FUNCTION AND CONSTRAINT VALUES 

C WRITE(7,SAVE) 

40 WRITE(7) IPRINT,NOV, ITMAX,NCDN,NSIDE, ICNDIR,NSCAL,NFDG, 
1F0CH,FDCHM,CT,CTMIN,CTL,CTLMIN,THETA,PHI,NAC,DELFUN,DABFUN, 
2LIN0BJ,ITRM,ITER, INFOG,IGOTO,INFO,OBJ, 

3RDUM, IDUM, 

4X,DF,6,ISC,IC,A,S,G1,62,C,MS1,B, VLB,VJB,SCAL,RSAV, ISAV,NCOUNT 
5,N1,N2,N3,N4,N5,ALPHAX, ABOBJl 
C WRITE(9, LINKF) 

REWIND 9 

WRITE(9, 50)NDV, (X( I), I»1,NDV) 

50 FQRMAT(I5/(6E13.5) ) 

C WRITE CONTROL CARD STORED IN PROCFILE GONOGO. 

101 F0RMAT(»1») 

IF( IGOTO.NE.O) GO TO 301 

CALL SRCH$$( INTS(2), ' GONOGO INTS ( 6) , INTS ( 8) , INTS ( 0 ) , ICQDE ) 
WRITE(12,101) 

CALL SRCH$$(INTS(4), ' GONOGO INTS ( 6) , INTS ( 8 ), INTS ( 0 ) , I CODE ) 

C WRITE ON TAPES IF GRADIENTS ARE REQUIRED 

REWIND 13 
JITER » JITER+1 

WRITE (13,14) JITER,NDV, NCON,(ISC(I ),I«1,NC0N) 

WRITE(14,45) 0BJ,(X(I),I»1,NDV),CT,CTL,(G(I),I*1,NC0N) 

301 IF(INF0.NE.2) GO TO 303 

CALL SRCH$$(INTS(2), 'STARTX*,INTS(6),INTS(4),INTS(0),IC0DE) 
WRITE(8,30) NDV, (X(I),I«1,NDV) 

CALL SRCH$$( INTS(4), »STARTX*,INTS(6),INTS(4), INTS ( 0) , ICODE) 


303 CALL SRCH$$(INTS(2), 'INFO', INTS(4),INTS(4),INTS(0), ICODE) 
WRITE(8,304) INFO 

304 FORMAT(Il) 

CALL SRCH$$(INTS(4), « INFOS INTS (4), INTS (4), INTS (0), ICODE) 

STOP 

END 
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APPENDIX C.- 


LISTING OF FRONT PROCESSOR 


PROGRAM F PFGSK INPUT# OUTPUT^ TAPE?# T APES- I NPUT#TAP66»DUTPUTJ 
SPAR FRONT PROCESSOR READS DESIGN VARIABLES AND 
PRINTS THEM IN SPAR SECTION PROPERTY FORMAT 
FUSELAGE MADE OF ROD BEAM MEMBRANE ELEMENTS 
DIMENSION X(50) 

NAMELIST/LINKF/NDV#X 
NOV-NUMBER OF DESIGN VARIABLES 
X(NDV)»DESIGN VARIABLES 

OV*S ARE X( D-SECTIONAL AREA OF STRINGER RODS 
X{2)»N0NDIM£NSI0NAL AREA OF BEAM 
X<3)-THICKNESS OF MEMBRANE PANEL 
READ(7#10) B10#B20»T0 
10 F0RMAT(3F10.3) 

C REA0(5#LINKF) 

READ(5#15)NDV# (X( I)#I*1#NDV) 

15 F0RMAT(I5/(6E13.5) ) 

XIN1»1./X(1) 

XIN2-1./X(2) 

XIN3»1./X(3) 

C WRITE E23 ELEMENTS 
WRITE(6#200) 

200 FORMAT!* E23 SECTION PROPERTIES*) 

WRITE(6#201) XINl 

201 FORMAT!' 1 *#F8.3) 

C COMPUTE VALUES FOR DSY CARDS 
C WRITE E21 ELEMENTS 

WRITE(6#202) 

202 FORMAT!* E21 SECTION PROPERTIES*) 

AREA0»!2. + B10 + B20)*‘T0 
AREA-AREA0>«'XIN2 

SC ALE -SORT! AREA /A RE AO) 

Bl-BIO^SCALE 

B2-820«SCALE 

T-TO*SCALE 

EI1-B1*(B2+2.*T)**3/12.-(B1-T)*B2^*3/12. 

ALPHAl-O. 

C-!B2 + 2.*T)*'Bl**2?2.-B2*!81-T)*‘!Bl + T)/2. 

C-C/AREA 

El2-2.n*Bl«*3/12.+2,*T*‘Bl*'! Bl/2 .-C)**‘2 
l + B2*T«*3/12. + B2«T*‘!C-T/2. )* + 2 
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ALPHA2-0. 

F«(2.*Bl+B2)*T+*3/3. 

Fl = Oo 

Zl» ( (Bl*B2)**2)>T/(A.’«'EIl) + C-T/2. 

Z2«0. 

THETA=0. 

Q1 = 0. 

Qe-’O. 

Q3«0. 

m-~{Bi-c) 

Y12».5*B2+T 

Y21=C 

Y22».5*B2+T 

Y31»C 

Y32*-( .5’»'B2+T) 

V41=-{Bl-rj 
Y^»2«-{ .5>^B2 + T) 

J« 1 

WRITE 5 6^103) ,) # E 1 A L PHAl , E 1 2> AL PH A2# AR E A 
WRITE (6.^ 100 3) Fl»Zlf Z2?THETA 

103 FnRfiATC'>DSY’’j»I2j»!)F12eAi.5 7„i! ) 

1003 FORMATdX^ 5F.l3»r)) 

WRITE (6/. 104) 01/»Q2j Q3; YU# Y12,Y21 

104 F0RMATC1K#6F12,4-. '5JM 

WRITE « 6# 1004 5 Y22# f 31# Y32.» Y4 1# Y42 

1004 FORMAT? lKj?5F12«4) 

C WRITE E41 ELEMENTS 

HRITE(6#300) XIN3 

300 format?* Shell section properties* /• i*#f8.3) 

WRITE(6#301 ) 

301 FnRMAT?’S*) 

STOP 

END 
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APPENDIX D.- LISTING OF END PROCESSOR. 

PROGRAM EPPUSaNPUT»OUTPl'T>7APE8,TAPE5»INPUTf TAPE6»QUTPUT) 
DIMEMSION A(4500)iB(lAOO) 

DIMENSION A1(A00),B1(400),C1( AGO) 

DIMENSION G(400) 

NAMELIST/EPIN/E23AL»E21AL>E41AL»NSE23»NSE21»NSE41 
C NAMELIST/LINKE/OR J#G 

REAO(5).EPINI 

CALL 0AL(4»11^A(1)*0>IEA»KADR>IERR>NWGS»NE»LB>ITYPE> 

14H0BJF>3HAUS»1» 1) 

OBJ*A(l) 

CALL 0AL{4>ll»A(l),0»IE4,KA0P>IERR»N*^DSfNE»LB#ITYPE^ 
14HSTRS/3HE23, 1 ) 

11 = 1 

INL-6+NSE23 
DO 2 IN»6»INL^6 
1=11 

Al( II ) «A( IN ) 

G{I)«A8S(A1(I1) )/E23AL-l. 

11=11+1 
2 CONTINUE 

CALL DAL(4fll^A(l),0»IEA,!<ADR»IERR,NWDS^NE,LB#ITYPE» 
14HSTRS*3HE21»1>1) 

Jl»l 

KNL1=(NSE21-1)’«‘52 + 6 
DO 4 IM»5#KNL1>52 
B1(J1)«ABS(A(IM)) 

A1(J1)=ABS(A(IM+1) ) 

I-Il+Jl-1 
GNUM»B1( Jl) 

IF ( Al( Jl) .6T.B1 ( Jl) ) GNUM.AKvII ) 

G(I )»GNUM/621 AL-1 . 

J1=J1+1 
4 CONTINUE 

CALL 0AL(4,11,A(1),0,IRA,KADR>IERR»NWDS^NE^LB»ITYPE, 
14HSTRS»3HE41, 1, 1) 

Kl = l 

LNLl»(NSE41+23) 

DO 10 NP=21,LNL1»23 
CKKl )=A(NP) 

B1 (K1 )»A (NP + 1 ) 

A1 (K1 )»A (NP + 2 ) 

1=11+ Jl+Kl-2 

GNUM«SQRT(Cl(Kl)* + 2 + Bl(Kl)+ + 2-Cl(Kl) + Bl(Kl)+3.*Al(Kl) + ’f2> 
G( I )»GNUM/£41AL-1 , 

K1=K1+1 
10 CONTINUE 

CALL FIN(OiO) 

REWIND 6 

C WRITE (6,LINKE ) 

WRITE(6#232) OBJ^G 

232 FORMATC^BEGIN C NMN I D« / E 1 3 . 5 / ( 6 El 3 . 5 ) ) 

WRITE(6^233) 

233 FORMAT(*END CNMNIO+) 

STOP 

END 
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APPENDIX E.- LISTING OF BLDELDS. 


PROGRAM BLDELDS (TAPE 5# TAPE2 3#TAP E 20> TAPE 21> TA PE 22» OUTPUT » 

THIS PROGRAM CREATES A RUNSTREAM THAT WILL CREATE 

DMDV AND DKDV FOR PARTICULAR ELEMENTS USED DESIGN VARIABLES 

DIMENSION EL(999),NSECT(999),TNAME1(8)#TNAME2(8),F0R(9) 

DIMENSION N0DVPE(999) 

DATA E21#E22»E23>EAl»EA3>EAA/3HE21>3HE22> 3HE23# 3HEA1,3HE43> 3HEA4/ 
DATA E31,E33/3HE31»3HE33/ 

DATA TNAMEl/AHDEF MHGD # AHGTIT, AHDIR ,AHNS »3+4HELTS/ 

DATA TNAME2/5+AH f AHNAME» AHNNOD>AHISCT/ 

DATA START# END# XNS ECT/ AH$ST A# AHSE ND# 4HNSEC / 

DATA YNSECT/3HNSE/ 

CALL SUBROUTINE TO REMOVE BEGINNING BLANKS 
CALL REMOVE 
READ INPUT 

NOEL-NUMBER OF ELEMENTS 

NODV-NUMBER OF DESIGN VARIABLE ELEMENTS 

VOR8»TYPE OF ANALYSIS (EX. BUCKLING) 

NDF»NUMBER OF DEGREES OF FREEDOM PER JOINT 
EL - ELEMENT NAMES CONTAINING DESIGN VARIABLES 
(EX, E21) 
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NSECT - LAST SECTION NUMBER USED FOR EACH DESIGN VARIABLE 
NODVPE - NUMBER OF DESIGN VARIABLES PER ELEMENT 

READ( 5,5) N0LC,N0DV, ISNOLC, J0INTS,NDF,N0EL,V0RB 

5 F0RMAT(6<1X,I^),1X,A4) 

RE AD (5,6) (EL(I ),NSECT(I),N0DVPE( I), I»1,N0EL) 

6 F0RMAT(6(1X,A3,1X,I3,1X, 13) ) 

NDF»NDF*NDF 

WRITE(20,2) 

2 FQRMAT(*CXQT TAB*/+ UPDATE»1*) 

LOOP ON NUMBER OF ELEMENTS 
DO 30 I * 1,N0EL 

CALL SUBROUTINE TO UPDATE TAB BY SETTING DESIGN VARIABLES 
TO UNITY. 

KCNT RETURNS THE NUMBER OF POSSIBLE DESIGN VARIABLES FOR 
A PARTICULAR ELEMENT. 

KCNT«99 MEANS THE ELEMENT NAME IS BAD 

CALL TABNPUT(EL(I),NSECT(I),KCNT) 

IF(KCNT.NE.99) GO TO 30 
PRINT 29,EL(I) 

29 FORMAT(* ELEMENT NAME ♦,A3,* DOES NOT EXIST*) 

GO TO 190 

30 CONTINUE 

CONCLUDE UPDATE AND SET UP DISABLES 

WRITE (20,31) 

31 F0RMAT(* UPDATE-0*) 

WRITE(20,32 ) 

32 F0RMAT(*CXQT DCU*/* COPY 1,2*) 

DO 35 I <= 1,3 

IF (I .GT.A) go to 33 
00 37 J a 1,N0EL 
TNAME2(I)-EL(J) 

WRITE (20, 34) TNAMEl ( I ) , TNAME2 ( I ) 

37 CONTINUE 
GO TO 35 

33 WRITE(20,34) TNAME 1 ( I ) , TNAME2 ( I ) 

34 FORMAT(* DISABLE 1, *, A4, IX, A4 ) 

35 CONTINUE 

SET COUNTER FOR TOTAL NUMBER OF DESIGN VARIABLES 

lELCNT = 0 
ICNTDV a 1 
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ISW«0 

ISAVCNT « 0 

READ IN RUNSTREAM AND CHECK FOR START OF A DESIGN VARIABLE 


39 

JCNT ■ 

0 





40 

READ(21#50) 

(FOR( J), 

J* 

1>9) 

50 

FORMAT 

(A4 

f A6 

#7A10) 




IF(EOF 

(21 

) ) 

170#60 



60 

IF(FOR 

(1) 

.EQ 

.START) 

GO 

TO 70 


WRITE( 

20> 

50) 

(FQR( J) 

r J 

«1>9) 


GO TO 

40 





70 

IFdSW 

.EQ 

.1 ) 

WRITE(20# 

75) 

75 

FORMAT 

(*C 

XQT 

ELD*) 




lSH-1 







READ(2 

1»80) 

(FOR( J)> 

J « 

1#9) 

80 

FORMAT 

(A3 

^A7 

>7A10) 




SET ICNT » NUMBER OF POSSIBLE DESIGN VARIABLES FOR AN ELEMENT 
ICNT«99 

DETERMINE POSSIBLE NUMBER OF DESIGN VARIABLES PER ELEMENT 

IF(F0R(1).EQ.E23.0R.F0R(1),EQ.E41.0R.F0R(1 )*E0.E4A) ICNT»1 
IF(FQR<1>,EQ.E31) ICNT«1 
IF(F0R(1),EQ.E21) ICNT«4 

IF( FOR(l) ,E0,E43.0R.F0R(1) . EO.E33) ICNT-12 
IF(F0R(1).EQ.E22) ICNT«21 
TNAME»FOR(l) 

IF(ICNT.EQ.99) TNAME»SAVNAM 
JCNT » JCNT+1 

IF(ICNT.NE.99) lELCNT-I ELCNT+1 

SET UNIT = 20 IF ONLY ONE DESIGN VARIABLE PER ELEMENT 
OTHERWISE SET UNIT » 22 (SCRATCH UNIT) 

IUNIT-22 
REWIND 22 

IF(ICNT.EQ.l) IUNir»20 

IF(ICNT.EQ.99.AND.ISAVCNT.EQ.l) IUNIT«20 

CHECK FOR REPEAT OF ELEMENT NAME 

IF(ICNT.NE.99) GO TO 86 
WRITE(IUNIT»85) SAVNAM 
85 F0RMAT(A3»77X) 

GO TO 860 
C 
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C READ DATA FROM UNIT 21 AND WRITE DATA ON UNIT 20 OR 22 
C DEPENDING UPON VALUE OF ICNT 
C 

86 WRITE (lUNIT, 50) ( FOR ( J ) , J -1 , 9 ) 

SAVNAM»F0R{1) 

ISAVCNT-ICNT 

60 TO 87 

860 ICNT » ISAVCNT 
GO TO 870 

87 READ(21^50) ( FOR ( J ) , 9) 

870 IF (FOR(l) ,EQ.XNSECT.0R.FDR( 1) .EQ. YNSECT) 60 TO 88 
IF(F0R(1 ) .EQ.ENO) GO TO 110 
WRITE(IUNIT,50) ( FOR ( J ) , J»1 , 9) 

GO TO 87 

88 IP1«NSECT(IELCNT)+JCNT 
WRITE(IUNIT#90) IPl 

90 F0RMAT(*NSECT»*>I3»71X) 

100 READ(21#50) ( FOR ( J ) ^ J *1> 9) 

IF<F0R(1).EQ.END) GO TO 110 
WRITE(IUNIT,50) ( FOR ( J ) , J»l, 9) 

GO TO 100 

SKIP THIS IF MORE THAN ONE DESIGN VARIABLE PER ELEMENT 

110 IF(ICNT.NE.l) 60 TO 120 

CALL SUBROUTINE TO CREATE REMAINDER OF RUNSTREAM 

CALL CRRS(ICNTDV#0#0>TNAME>NDF>NODVPE(IELCNT) ) 

ICNTDV » ICNTDV+NODVPE( lELCNT) 

GO TO 39 

LOOP ON POSSIBLE NUMBER OF DESIGN VARIABLES PER ELEMENT 
READ FROM UNIT 22 

WRITE ON UNIT 20 

120 DO 160 I ■ 1, ICNT 

IF(I.NE.l) WRITE{20,75) 

REWIND 22 

130 READ(22#50) ( FOR ( J ) # J»1 » 9) 

IF(E0F(22)) 150#1A0 

lAO IF(F0R(1) .EQ.XNSECT.OR.FORI 1 ).EQ. YNSECT) GO TO 141 
WRITE(20#50) (F0R( J)»J«1#9) 

60 TO 130 

141 WRITE(20#90) IPl 

GO TO 130 
C 

C CALL SUBROUTINE TO CREATE REMAINDER OF RUNSTREAM 
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C 

150 CALL CRRStICNTDV# ICNT# I>TNAME>N0F#N0DVPE( lELCNT) ) 

IP1»IP1+1 
160 CONTINUE 

ICNTOV » ICNTOV+NOOVPE( lELCNT) 

GO TO 39 

170 WRITE(20#180) 

180 FORMAT!* TOC 2*/*CXQT EXIT*) 

190 STOP 
END 

SUBROUTINE TABNPUT ( EL N AME» NSCT, I CNT ) 

C 

C THIS SUBROUTINE CREATES TAB PROCESSOR INPUT FOR A RUNSTREAM 
C DEPENDING UPON ELEMENTS USED. 

C ALL DESIGN VARIABLES ARE SET TO UNITY. 

C ICNT » NUMBER OF DESIGN VARIABLES FOR A PARTICULAR ELEMENT 
C ICNT » 99 MEANS THE ELEMENT NAME IS BAD 
C 

DIMENSION ELEMENT(21) 

DATA E21>E22#E23>E41#E43#E4A/3HE21#3HE22>3HE23>3HEA1,3HE43>3HE44/ 
DATA E31>E33/3HE31»3HE33/ 

DATA BA>BB>BC»SA»SB/2HBA^2HBB#2HBC>2HSA»2HSB/ 

IF(ELNAME.EQ.E21) ELID=BA 
IF(ELNAME.EQ.E22) ELID=8B 
IF(ELNAME.EQ.E23) ELID»8C 

IF(ELNAM£.EQ.E41.0R.ELNAME.EQ.E43.0R.ELNAME.EQ.E3l.0R. 

1 ELNAME.EQ.E33) ELID=SA 
IF(ELNAME.EQ.E44) ELID-SB 
WRITE(20,10) ELID 
10 FDRMAT(2X#A2) 

ICNT=99 

IF(ELNAME.EQ.E23.0R.ELNAME.EQ.E41.0R.ELNAME.EQ.E44) GO TO 30 
IF(ELNAME.E0.E31 ) GO TO 30 
IF{£LNAME.EQ.E21) GO TO 50 

IF(ELNAME.EQ.E43.0R.ELNAME.EQ.E33) GO TO 100 
IF (ELNAME.EQ.E22) GO TO 150 
PRINT 20jELNAME 

20 FORMAT!* ELEMENT NAME *>A3#* DOES NOT EXIST*) 

GO TO 210 

ELEMENT NAMES E23 > E31 , E41 , E44 
ICNT = 1 

30 K«NSCT+1 

WRITE!20,40) K 
40 F0RMAT!1X#I3#* 1.0*) 

ICNT»1 
GO TO 210 
C 
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C ELEMENT NAME E21 



ICNT 

» A 

50 

00 90 I » 

1#^ 


DO 60 J = 

1#10 


ELEMENT! J ) 

« 0.0 

60 

CONTINUE 



IF(I.EQ.l) ELEMENT(1)»1.0 
IF(I.EQ,2) ELEMENT(3)»1.0 
IF(I.EQ.3) ELEMENT(5)«1.0 
IF(I.EQ.A) ELEMENT(6)»1.0 
K-I+NSCT 

WRITE { 20# 70) K# ( ELEMENT ( J )# J= 1# 10 ) 
70 F0RMAT(>> DSY ♦# II# 10( 1X#F3 .1 ) ) 
WRITE( 20#80) 

80 FnRMAT{2X#ll(<'0.0 ♦)) 

90 CONTINUE 
ICNT»A 
GO TO 210 
C 

C ELEMENT NAMES E33 # EA3 
C ICNT » 12 

C 

100 DO lAO I ■ 1#12 
DO 110 J » 1#12 
ELEMENT( J )»0.0 
110 CONTINUE 

ELEMENT(l)»l.E-06 

ELEMENT(3)»l.E-06 

ELEMENT(6)=1.E~06 

ELEMENT(7)=l.E-06 

ELEMENT(9)»l.E-06 

ELEMENT(12)»l.E-06 

ELEMENT(I)»1. 

K-I+NSCT 

IF(I.EQ.l) WRITE(20#115) 

115 FORMAT!* FORMAT ^UNCOUPLED* ) 

WRITE! 20#120) K 

120 F0RMAT!2X#I2#1X#9!*1.0 *)) 

WRITE!20#130) ! EL EMENT ! J ) # J »1# 6) 
130 F0RMAT!1X#6!1X#E10.2) ) 

WRITE! 20# 130) ! E LEME NT ! J ) # J *7# 12 ) 
lAO CONTINUE 
ICNT»12 
GO TO 210 
C 

C ELEMENT NAME E22 
C ICNT - 21 

C 
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150 00 200 I » 1,21 
00 160 J • 1,21 
ELEMENTU)*0.0 
160 CONTINUE 
K»I+NSCT 
ELEMENTCI)«1. 

MRITE<20,170) K,ELEMENT(1) 

170 F0RMAT(2X,I2,1X,E10.2) 

KK-2 

DO 190 J ■ 2,6 
L»KK+J-1 

WRITE(20,180) (ELENENT(N),N»KK,L ) 

180 F0RMAT(2X,6(F3.1,1X) ) 

KK-KK+J 
190 CONTINUE 
200 CONTINUE 
ICNT-21 
210 RETURN 
END 

SUBROUTINE CRRSCICNTDV, ICNT,J,TNAME,NOF,NOOVEl) 

THIS SUBROUTINE CREATES REMAINDER OF RUNSTREAM 
TO FIND DMDV AND DKDV. 

USE NAMES 

OMDV DIAG 0 I AND 
OKOV SPAR 25 I 

WHERE I « 1 TO NUMBER OF POSSIBLE DESIGN VARIABLES 

DIMENSION BA(A), 88(21), SAEL ( 12 ), TNAMEl { 9 ) , TN AME2 ( 9 ) 

DATA DK,OM,CK,CM,SK» SM/2H0K, 2HDM, 2HCK,2HCM,2HSK, 2HSM/ 
DATA DV/2HDV/ 

DATA BB/2H11,2H21,2H2 2,2H31,2H32,2H33,2H41,2H-(»2,2HA3, 

1 2HAA,2H51,2H52,2H53,2H5A,2H55,2H61,2H62, 2H63, 2H64, 2H6 5, 

2 2H66/ 

DATA SAEL/2H11,2H12,2H2 2,2H13,2H23,2H33,2H4<»,2H45, 2H55, 

1 2HA6,2H56,2H66/ 

DATA OEM/IOHOEM DIAG 0/ 

DATA KSP/7HK SPAR ( 

DATA CHA/9HCHANGE 2,/ 

DATA C0P/9HC0PY 1,2 / 

DATA BA/2HIX,2HIV,2HDA,2HJ0/ 

DATA TNAME1/4H0EF ,4HGD , 4HGT IT, 4HELTS, 4HNS ,4HKMAP, 

1 4HAMAP,4HMASK,4H0IR / 

DATA TNAME2/3+4H ,4HMASK,3’»'4H ,4HEFIL,4H / 

DATA DM0V,DKDV,SPAR/4HDMDV,4H0KDV,6H SPAR / 

DATA DIAG/6H DIAG / 

JCNTOV • ICNTDV 
C 

C XQT E , EKS , TOPO , K , DCU 
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c 

WRITE(20»10) 

10 F0RMAT(*CXQT E*/*CXQT EKS+/>4<CXQT T0P0*/*CXQT K*/*CXQT DCU*) 

DISABLE DATA SETS 

DO 30 I « 1>9 
IF(I,EQ.8) TNAMEK8) »TNAME 

IFn.EO,1.0R,I,EQ.2.0R. I.EQ.3.0R .1 .EQ.9) TNAME2 ( I ) =TNAME 
HRITE(20,20) IN AMEU I ) , TNAME2 ( I ) 

20 F0RMAT(« DISABLE 1»*A4>2X,AA) 

30 CONTINUE 

SET UP ELEMENT NAMES FOR CHANGE AND COPY STATEMENTS 

DDVDV 
DDK-DK 
DDM-DM 

IF(ICNT.NE.A) GO TO 31 
DDV»BA( J) 

GO TO 39 

31 IF(ICNT.NE.12) GO TO 32 
00K«CK 
DDM-CM 
DDV«SAEL(J) 

GO TO 39 

32 IF(ICNT.NE.21) 60 TO 39 
DDK»SK 
DDM»SM 
DDV-BB(J) 

39 DO 80 UK » 1,N0DVEL 
WRITE(20,40) CQP>DEM 
AO F0RMAT( A9,A10,>«' 0*) 

MRITE(20j50) CHA>0EM» D0M>D0V>DIAG# JCNTDV 
50 F0RMAT(A9^ AlO# ♦ 0» ♦» 2A2» A6> *0 *#I3) 

WRITE(20,60) COPfKSPjNDF 
60 FQRMAT( A9,A7> I2>>»' 0*) 

WRITE (20, 70) CHA,KSP,NDF,OOK,DD\/,SPAR,NDF, JCNTDV 
70 F0RMAT(A9,A7,I2,+ 0, ♦, 2A2, A6, 1 2, I A ) 

JCNTOV » JCNTDV+1 
80 CONTINUE 
RETURN 
END 

SUBROUTINE REMOVE 

THIS SUBROUTINE REMOVES THE LEADING BLANKS FROM EACH 
LINE IN THE RUNSTREAM. 

DIMENSION 0ATIN(80) 


72 



DATA BLAnK/IH / 

REWIND 23 

1 READ(23#2I ( DAT IN( I ) , I / 80 ) 

2 F0RMAT(80A1) 

IF(E0F(23)) 7,3 

3 IF(DATIN<1).NE. BLANK) GO TO 50 
DO 6 I » 2,80 

IF(DATIN< I) .NE.BLANK ) GO TO 60 

6 CONTINUE 

60 L«I-1 

K-80-I 

00 <» J » 1,K 
OATIN<J) « DATIN(J*L» 

A CONTINUE 
K»K + 1 

DO 5 J « K,80 
OATIN(J) » BLANK 
5 CONTINUE 

50 WRITE(21,2) ( DATIN ( J ) , J » 1, 80 ) 

GO TO 1 

7 REWIND 21 
RETURN 
END 


73 



o o o o 


APPENDIX F.- LISTING OF GRDXXXX 


PROGRAM GN6RDRS(lNPUT#TAPE3 0»TAPE31>TAPei0, TAPE 5* INPUT, OUTPUT) 
C 

C THIS PROGRAM CREATES A REPEATABLE SPAR RUNSTREAM 
C FOR CALCULATING DERIVATIVES WITH RESPECT TO 
C DESIGN VARIABLES. THE SIZE OF THE RUNSTREAM VARIES 
C WITH THE NUMBER OF LOAD CASES (NOLO AND THE 
C NUMBER OF DESIGN VARI ABL ES ( NOD V ) . 

C 

C THE RUNSTREAM IS OUTPUT ON UNIT 10 
C 

DIMENSION EL(999),NSECT(999 ),NUM<8),MFQRM(5),N0DVPE(999) 

DATA NUM/1H1,1H2,1H3,1H4,1H5,1H6, 1H7,1H8/ 

DATA BUCK,VIBR/AHBUCK, ^HVIBR/ 

DATA E21,E22,E33,E<t3/3HE21,3HE22, 3HE33,3H£A3/ 

DATA MF0RM(1)/10H(*Z11=UNI0/ 

DATA MF0RM( A)/10Hm»lX,*Z^/ 

DATA MFQRM(5)/9H),I1,>5'H')/ 

DATA MF0RM7/9HN(Z10,Z«, / 

DATA MF0RM8/5HN(Z*,/ 

READ INPUT 

NOEL ■ NUMBER OF DIFFERENT ELEMENTS 
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NOLC » NUMBER OF LOAD CASES 

NODV « TOTALNUMBER OF DESIGN VARIABLES 

ISNOLC - STARTING NUMBER FOR DERIVATIVE LOAD CASES 

NDF • NUMBER OF DEGREES OF FREEDOM PER JOINT 

NORB - TYPE OF ANALYSIS (EX. BUCKLING) 

JOINTS » NUMBER OF JOINTS IN THE MODEL 
NODVPE » NUMBER OF DESIGN VARIABLES PER ELEMENT 
EL a NAMES OF ELEMENTS CONTAINING DESIGN VARIABLES 
(EX. E21) 

READ(5>60) NOLC » N00V» ISNOLC > JOINTS# N0F» NOE L>V0RB 

60 F0RMAT(6(1X#I4),1X#A4) 

READ(5#61) (EL(I)#NSECT(I)#N0DVPE(I)#I*1#N0EL) 

61 F0RMAT(6(1X,A3#1X»I3»IX#I3) ) 

ISN0LC«ISN0LC~1 

NOF»NDF«NDF 

DETERMINE IF THERE IS A E21 , E22 # E33 » OR EA3 ELEMENT 

ICNTDV « 1 
00 168 I « 1#N0EL 
JJ ■ NODVPE (I) 

DO 160 J ■ 1#JJ 
IF(EL(I).EQ.E21) GO TO 161 
IF(EL(I) .EQ.E22) GO TO 163 
IF(EL(I),EQ.EA3) GO TO 16A 
IF(EL(I).EQ.E33) GO TO 165 
GO TO 166 

CALL SUBROUTINE TO FIND DK/DV AND DM/DV FOR E21 ELEMENTS 

161 CALL 0KDVE21(ICNTDV#HDF) 

GO TO 166 

CALL SUBROUTINE TO FIND DK/DV AND DM/DV FOR E22 ELEMENTS 

163 CALL DKDVE22( ICNTDV#NDF) 

GO TO 166 

CALL SUBROUTINE TO FIND DK/DV AND DM/DV FOR EA3 ELEMENTS 

16<r CALL DKDVEA3( ICNTDV#NDF ) 

GO TO 166 

CALL SUBROUTINE TO FIND DK/DV FOR E33 ELEMENTS 

165 CALL DKDVE33(ICNTDV#NDF) 

166 ICNTDV ■ ICNTDV+1 
160 CONTINUE 
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168 


CONTINUE 

IF(VORB,EQ.VIBR) GO TO 320 
WRITE(10>172) 

172 F0RMAT(*CXQT AUS+) 

FIND OBJECTIVE FUNCTIONS 
162 WRITE(10,175 ) JOINTS 

175 F0RMAT(+ SYSVEC;UNIT MEC*/* I»lJ J»1,*,I8>«; 1,0*/ 
1 ♦ DEFINE UN-UNIT VEC*) 

DO 200 I « IjNODV 
WRITEtlO^lBO) I, I 

180 F0RMAT(* DEFINE W*#I3#*»DMDV DIAG 0 *>I3) 

200 CONTINUE 

00 250 I » 1,N0DV 
WRITE(10^220) I, I 

220 F0RMAT(* OBJF G*#I3#* 1 1»XTY(UN> W*>13>*) ♦ ) 

250 CONTINUE 

WRITE(10^260) 

260 FORMAT(*CXQT DCU*) 

DO 300 I « 1»N0DV 
HRlTE(10f270) I 

270 F0RMAT(* PRINT 1 OBJF 6*>I3» 

300 CONTINUE 

CALL SUBROUTINE TO CREATE RUNSTREAN FOR DERIVATIVE 

CALCULATION 


FIND APPLIED FORCES AND MOMENTS 
WRITE<10,2) 

2 FORMAT(*CXQT AUS*/* OUTLlB-3*) 

DO 4 I » 1,N0LC 
WRITE(10#3) I, I 

3 FORMAK* DEFINE F*»I2>*-STAT DISP *>I2* 1*) 

4 CONTINUE 

DO 6 I » 1,N0DV 
WRITE(10>5) I#NDF,I 

5 FORMAT!* DEFINE L*>I3»*-DKDV SPAR *»I2#I4) 

6 CONTINUE 

DO 10 I ■ 1,N0LC 
NWNOLC-I+ISNOLC 
MRITE(10>7) NWNOLC 

7 FORMAT!* ALPHA; CASE TITLE *#I3) 

DO 9 J » 1>N0DV 

WRITE!10,8) J#I#J 

8 F0RMAT!1X#I3»* "LOAD CASE *,I2#* DERIVATIVE 1 DESIGN VARIABLE *, 
1 13) 
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9 CONTINUE 
10 CONTINUE 

IFINQOV.NE.ll WRITEdO, 2003 I 
2003 FORMAT!* OUTLIB-^*) 

IFINOOV.EQ.l) WRITE(10#6006) 

6006 FORMAT!* QUTLIB»3 ♦) 

DO 1010 I « 1>NQLC 

NWNOLC»I+ISNOLC 

DO 1005 J « 1»N0DV 

WRITE !10d003) NWNQLC#J#J»I 

1003 FORMAT!* APPL FORC *»I3,I3,*» PR00UCT(-1.0 L*> 1 3# 1 • 0 F*> I 2# * ) * ) 
1005 CONTINUE 
1010 CONTINUE 

DO 15 I ■ 1»N0LC 
NWNOLC-I+ISNQLC 
IF!N0DV.E0.1) GO TO 50 
IF!I,EQ.l) WRITE!10,2009) 

2009 FORMAT!* INLIB»<f*/* 0UTLIB«3*) 

DO 13 J ■ 1#N0DV>9 
IT0P»9 

ICHK»N0DV-J+1 
IF!ICHK.IT.9) ITOP-ICHK 
DO AO K » 1#IT0P 
L»K+J-1 

WRITE!10»12) K#NWNOLC#L 

12 FORMAT!* DEFINE Z*»I1»*=APPL FORC *,I3»IX,I3) 

AO CONTINUE 

IF!ITOP.EQ.l) GO TO A1 
MF0RM!2)»MF0RM7 
IFU.EQ.l) MF0RM(2)»MF0RH8 
MF0RM!3)«NUM!ITDP-1) 

WRITE!10»MFQRM) !K»K=1>ITQP) 

A1 IF !ITOP.EQ.l) WRITEdO, A2> 

A2 F0RMAT!*Z11=UNI0N!Z10 Zl)*) 

MRITE!10, A5) 

A5 FORMAT!* INLIB=3*/*Z12»UNI0N! Zll ) */♦ DEFINE Z10-Z12*/ 

1 * INLIB«A*) 

13 CONTINUE 
WRITE!10»A6) NWNQLC 

A6 FORMAT!* APPL FORC *,I3,* 1 «UN ION ! ZIO )* ) 

15 CONTINUE 

FIND STRESS AND DISPLACEMENT DERIVATIVES 

50 DO 32 I * 1,N0LC 
WRITE!10,16) NODV 

16 FORMAT!*CXQT SSOL*/* RESET L 1-1, L2»*, I 3) 

NWNOLC-I+ISNOLC 

WRITE!10,17) NWNQLC 
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17 

FORMAT!* RESET QLIB»3*/* RESET SET**, 13) 
WRITE(10,18) 



18 

FORMAT(*tXQT VPRT*/* LIB«3*) 
WRITE(10,19) NWNOLC 



19 

FORMAT(* PRINT APPL FORC *,I3) 
WRITE(10,20) NWNOLC 



20 

FORMAT!* PRINT STAT DISP *,I3) 
IF(VORB.EQ.BUCK) GO TO 400 
WRITE!10,21) NODV 



21 

FORMAT!*CXQT GSF*/* RESET L 1* 1, L2 **, 1 3/ * 
WPITE!10,22) NWNOLC 

RESET 

QLIB*3*) 

22 

FORMAT!* RESET SET**, 13) 
WRITE(10,23) NODV 



23 

FORMAT!*CXQT PSF*/* RESET Ll*l, L2-*, 13/* 
WRITE!10,24) NWNOLC 

RESET 

QLIB»3* ) 

24 

FORMAT!* RESET SET**, 13) 
GO TO 25 




SET UP RUNSTREAM FOR BUCKLING ANALYSIS 

400 IDV = 0 

DO 470 J » 1>NQEL 
NODVEL « NODVPE(J) 

DO 460 K » 1>N0DVEL 
IDV » IDV+1 
WRITE(10#425) 

425 FORMAT(»CXQT GSF*/+ RESET EMBEDal*) 

WRITE(10#430 ) IDV>IDV#NWNOLC 

430 FORMATC* RESET L 1 = *# I 3> L 2 »*» 1 3> ♦ > SET»* , I 3, ♦ » QL I B » 3* ) 
HRITEdOf 435) 

435 FORMATC^tXQT PSF*) 

HRITE(10#430) IDV> IDV#NyNOLC 
WRITE(10>440) 

440 FORMAT(*CXQT KG+) 

WRITE(10>445) 

445 FORMAT(<'CXQT DCU*) 

WRITE(10>450) NDF»NDF,IDV 

450 FORMATC* CHANGE 1>KG SPAR ♦>I3#* 0>DKG SPAR ♦>13,14) 
WRITE(10,455) NDF,IOV 
455 FORMAT(* COPY 1,3 DKG SPAR *,13, 14) 

460 CONTINUE 
470 CONTINUE 

25 WRITE (10,251) 

251 FORMAT(<‘CXQT DCU*) 

WRITE(10,26) NWNQLC,NWNOLC 

26 FORMAT(* CHANGE 3,STAT DISP ♦,I3,* 1,DDIS DISP *,I3,« !♦ ) 
CHANGE DATA SET NAMES 
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IDV»0 

00 31 J-1»N0EL 
NOSECT « NSECT(J) 

NODVEL « NODVPE(J) 

DO 30 K«1^N0DVEL 
DO 271 U « 1#N0SECT 
IDV»IDV+1 
DO 27 KK • 1,N0EL 
IF(EUKK) .EQ.E21) GO TO 29 
IF(EL (KK) .EQ.E22) GO TO 29 
IF(EL(KK) .EQ.E43J GO TO 29 
IF( EL<KK) .EQ.E33) GO TO 29 

WRITE <10 >28) EL(KK)>NWNOLC> IDV> EL ( KK ) > NWNQLC# IDV 
28 FORMAT(* CHANGE 3>STRS * > A3 > 21 DSTR ♦»A3>2I4) 
HRITE(10>285) EL ( KK ) > NWMOLC > IDV 
285 FORMAT!* COPY 3>4 DSTR *>A3>2I4) 

GO TO 27 

STORE BEAM CROSS DERIVATIVES 


29 WRITE(10>290) E L < KK ) > NWNOLC > IDV, E L ( KK ) > NWNOLC » IDV 
290 FORMAT!* CHANGE 3>STRS *> A3 > 2 1 4> *> DF AM *>A3>2I4) 

WRITE!10>292) EL ! KK ) > NWNQLC , I DV 
292 FORMAT!* COPY 3>4 OFAM *>A3>2I4) 

27 CONTINUE 
271 CONTINUE 

30 CONTINUE 

31 CONTINUE 

32 CONTINUE 

320 IF!VORB.EQ.VIBR) WRITE ! 10> 260 ) 

DO 333 I » 1»N0DV 
WRITE! 10> 332) I 

332 FORMAT!* COPY 1>4 OBJF G*>13) 

333 CONTINUE 
WRITE !10>33) 

33 FORMAT!* TOC 3*/* TOC 1*/* TOC 4*/*tXQT EXIT*) 
STOP 

END 
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APPENDIX G.- LISTING OF DRVXXXX 


PROGRAM 0RVSTRS(INPUT«65j»TAPE30#TAPE31#TAPE5»INPUT 
1 ,TAPE15,TAPE16,0UTPUT=65,TAPE6=0UTPUT) 

THIS PROGRAM COMPUTES THE STRESSES AND STRESS DERIVATIVES 
FOR SPAR BEAM ELEMENTS USING SPAR LIBRARIES AS INPUT. 

THE STRESSES ARE COMPUTED FROM SPARLA USING DATA SET 
FAMS MASK I 1 WHERE I = 1 TO NOLC 
THE STRESSES ARE WRITTEN BACK ONTO SPARLA USING DATA SET 
STRS MASK I X WHERE I « X TO NOLC 

THE STRESS DERIVATIVES ARE COMPUTED FROM SPARLC USING DATA SET 

DFAM MASK I J WHERE I » ISTRTLC TO ISTRTLC+NOLC AND J * X TO NODV 
THE STRESS DERIVATIVES ARE WRITTEN BACK ONTO SPARLC USING DATA SET 

DSTR MASK I J WHERE I » ISTRTLC TO ISTRTLC+NOLC AND J = 1 TO NODV 

COMMON KORE>KEVEN,Z(l) 

DIMENSION NAMEK2) »NAME2(2)>FL (999),N0DVPE(999) 

DATA E2I#E22f E33» E43/3HE21>3HE22# 3HE33#3HE43/ 

DATA NAME1/^HFAMS>4HDFAM/ 

DATA NAME2/4HSTRS>^HDSTR/ 

READ INPUT VALUES 

NOEL IS THE NUMBER OF DIFFERENT ELEMENTS 
NDF IS THE DEGREES OF FREEDOM PER JOINT 
VORB IS THE TYPE OF ANALYSIS (EX. BUCKLING) 

JOINTS IS THE NUMBER OF JOINTS IN THE MODEL 
NOLC IS THE NUMBER OF LOAD CASES 
NODV IS THE NUMBER OF DESIGN VRIA3LES 

ISTRTLC IS THE STARTING NUMBER FOR THE DERIVATIVE LOAD CASES 
EL *NAMES OF ELEMENTS USED AS DESIGN VARIABLES 
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C CEX. E21) 

C 

READ(5fl) NOLC^NODV, I SNOLC # JO INTS »N0F #N0EL # V0R8 

1 F0RMAT<6aX,IA)»IX,A4) 

READ(5>2» (EL(I)#NSECT#N0DVPE(I),I*1,N0EU 

2 F0RMAT(6(1X# A3^ 1X,I3#1X> I3>) 

ISTRTLC«ISN0LC"1 
CALL RSET(U>0#0) 

LOOP ON THE NUMBER OF LOAD CASES 

DO A NCNTLC » 1>NQLC 

LOOP ON NUMBER OF ELEMENTS 

DO 10 ISTRS = 1>N0EL 

CHECK FOR E21#E22#E33#EA3 ELEMENT 

IF(EL (ISTRSJ .EQ.E21,DR.EL(ISTRS).EQ.EA3.0R. 

1 EL( ISTRS). EQ.E22. OR. EL( ISTRS). EQ.E33) GO TO 9 
GO TO 10 

9 NODVEL » NOOVPE( ISTRS) 

SET NUMBER OF STRESS OR STRESS DERIVATIVES TO BE WRITTEN 

IF(EL(ISTRS).EQ.E21.0R.EL(ISTRS) .EQ.E22) IWRDCNT « 8 
IF(EL(ISTRS).EQ.EA3.0R.EUISTRS).EQ.E33) IWRDCNT » A 

LOOP ON NUMBER OF DESIGN VARIABLES PER ELEMENT 

DO 90 JK s 1, NODVEL 
REWIND 15 
REWIND 16 

READ FAMS MASK FROM SPARLA AND COMPUTE STRS MASK 

CALL RDATSET(A, NAME1(1)» NCNTLC >1,NE1,1»0# EL ( ISTRS) ) 
REWIND 15 
REWIND 16 

CALL SUBROUTINE TO WRITE STRESSES ON SPAR LIBRARY 


WRITE STRS MASK ON SPARLA 

CALL WRTDATA(A,NAME2(1)#NCNTLC»1#NE1»EL(ISTRS)>IWRDCNT) 
90 CONTINUE 
10 CONTINUE 
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ICNTOV « 1 


LOOP ON NUMBER OF ELEMENTS 
DO 3 ISTRS=1>N0EL 

IF (EL ( ISTRS ) .EQ.E21.0R.EL( ISTRS) .EQ.E43.0R. 

1 EL(ISTRS).EQ.E22.QR.EL(ISTRS).EQ.E33) GO TO 19 
GO TO 3 

19 DO 30 ICNTOV » 1»N0DV 
I8EAM * 0 

NODVEL » NODVPE ( ICNTOV) 

SET SWITCH FOR BEAM ELEMENT 

IF (EL ( ICNTOV) .EQ.E21 .OR .EL ( ICNTOV ).E0.E22) IBEAM-1 
N3»NCNTLC+ISTPTLC 

SET NUMBER OF STRESSES AND STRESS DERIVATIVES TO BE WRITTEN 

IF(£L(ICNTDV).EQ.E21) IWRDCNT * 8 
IF(EL( ICNTOV) .EQ.E43) IWRDCNT = A 

LOOP ON NUMBER OF DESIGN VARIABLES PER ELEMENT 

DO 20 JK » 1>N0DVEL 

READ DFAM MASK FROM SPARLC AND COMPUTE DSTR MASK 
REWIND 15 

CALL RDATSET(4,NAME1(2),N3^ ICNTOV^ NE1» 0> I BEAM> E L ( ISTRS )) 
REWIND 15 

WRITE DSTR MASK ON SPARLC 

CALL WRTDATA( 4, NAME2(2 ) >N3»ICNTDV# NE If EL ( I ST P S ) > I WR OCN T) 

20 CONTINUE 
30 CONTINUE 
'3 CONTINUE 

4 CONTINUE 

CALL FIN(OfO) 

STOP 

END 

SUBROUTINE WRTDATA(NUfNlfN3flSTRSfNElfX2» IWRDCNT) 

SUBROUTINE TO WRITE STRESS AND STRESS DERIVATIVE 
DATA SETS BACK INTO SPAR LIBRARIES 

COMMON KOREfKEVENfZ(l) 
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C SET UP BLOCK SIZES 

C IF NWD3 GT OPEN CORE SIZE > USE MORE THAN ONE BLOCK 
C 

NW03-NE1*IWRDCNT 
L83-NW03 

IF(NWD3.LT.K0RE) GO TO 61 
LB3*NW03/2 

61 ISW-0 

IF(NW03.LT.K0RE) ISW-1 

CALL 0AL(NU»0>Z(1)>K0RE/1»KA0R3,IERR,NWD3>NE1>LB3»-1> 

1 N1>X2#N3,ISTRS) 

62 Il-l 

63 I2»I1+IWRDCNT-1 
IF(I2.GT4LB3) go to 66 

READ STRESSES OR STRESS DERIVATIVES OFF UNIT 15 

READUS) (Z(I ), I-Il^ 12) 

IF(E0F(15)) 66>65 

65 Il-Il-MHRDCNT 
GO TO 63 

WRITE STRESSES OR STRESS DERIVATIVES ONTO SPAR LIBRARY 

66 CALL RI0(NU»10,2#KADR3/Z(1) ^LS3) 

IF(ISW.EQ.l) GO TO 67 
ISW»1 

LB3-NW03-LB3 
GO TO 62 

67 RETURN 

end 

SUBROUTINE RDATSET(NL',N1,N3, ISTRS,NE1> IBEAM> X2) 

SUBROUTINE TO READ SPAR L I 3RARY, STOP E DATA, AND 
COMPUTE STRESSES OR STRESS DERIVATIVES 

COMMON KORE^KEVEN,Z(l) 

DATA E21, E22> E33/E43/3HE21, 3H6 22,3HE3 3^3HEA3/ 

SET UP BLOCK SIZE 

CALL DAL(NU^10> Z( 1) »KORE>l>KA0Ri>IERR»NV01»NEl»L31> ITYPE>N1> 
1 X2>N3>ISTPS) 

NIlaLBl/NEl 
KCNT»0 

NBLK1«NW01/LB1 

IF (N9LK1+L31 ,NE .NWOl) N SLK 1 =N B L'< 1 + 1 
00 6 J = 1>MBLK1 
NL31«LB1 

IF(J.EQ«NBLK1) NL31»NW01-(N9LKl-l)*L31 

READ DATA FROM SPAR LIBRARY 

CALL Ria(NU»20^2^KA0Rl»Z(l)^NL31> 

DO 30 JCNT»1» NLBlrNIl 
KCNT«KCNT*1 

CALL SUBROUTINE TO COMPUTE STRESSES AND STRESS DERIVATIVES 

IF(X2.EQ.E21.0R .X2.EQ.E22) CALL BM STR S ( I S W # KC NT» JCN T# I BE AM ) 
IF (X2.EQ.E^3.0R.X2.EQ,E33) CALL P L TSTR S ( I SW> KCNT> JCNT ) 

30 CONTINUE 
6 CONTINUE 
RETURN 

END 83 



APPENDIX H.- LISTING OF DKDVE21. 


SUBROUTINE DKDVE21(NDVJIN»NDF ) 

DIMENSION X(20) 

NAMELIST/LINKF/NDViX 

C 

C THIS SUBROUTINE CREATES A SPAR RUNSTREAM TO CALCULATE 
C 


c 

DK 

OK 

DA 

DK 

DI DK 

DI 

DK 

DJ 

c 

c 



t ^ 

X 

Y 

0 

c 

DV 

DA 

DV 

DI 

DV DI 

DV 

DJ 

DV 

c 

I 


I 

X 

I Y 

I 

0 

I 

c 

c 

c 

c 

DM 

DM 

DA 

DM 

01 DM 

DI 

DM 

DJ 

c 

c 


■ 

. ♦ . 

X 

1 

1 

1 

4- 

>- 1 
1 
I 

0 

c 

DV 

DA 

DV 

DI 

OV DI 

DV 

OJ 

DV 

c 

I 


I 

X 

I Y 

I 

0 

I 


C 

WRITEdO,!) 

1 FORMAT(*CXQT AUS*) 

WRITE(10#3)NDFj NDVJIM>N0F,N0VJIM 

3 FORMAT(» DEFINE A1»DKDA SPAR ♦#I2>H/+ DEFINE A2-DKIX« 
1 ♦ SPAR ♦,I2,IA) 

WRITE(10,A)NDF,NDVJIM#NDF»NDVJIM 

4 FORMAT(* DEFINE A3=DKIY SPAR <‘>I2>I4/* DEFINE A4»DKJ0* 
1 ♦ SPAR ♦#I2,I4) 

VJRITE(10>203)ND\/JIM#NDVJIM 

203 FORMAT!* DEFINE B1=DMDA DIAG 0 *#13/* DEFINE B2»DMIX* 

1 ♦ DIAG 0 *^I3) 

WRITE (lOf 204)NDVJIM»NDVJIM 

204 FORMAT!* DEFINE 83 = DMIY DIAG 0 *d3/* DEFINE B4-DMJ0* 

1 * DIAG 0 *fl3) 

C 

C COMPUTE DA/DV 
C 

DADV*1 . 

C 

C READ IN AND SET UP INITIALIZATION VALUES 
C 

C READ IN CONSTANTS FROM UNIT 30 
C 

READ!30#5) B10#B20fTO 

5 FQRMAT!3F10.3) 

C 

C READ IN DESIGN VARIABLES FROM UNIT 31 
C 

READ! 31, LINKF ) 
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AREAO-CZ.fBlO + eaO^TO 
AREA»AREAO/X( 2) 

SCALE»SQRT(AREA/AREAO) 

Bl-BlO'fSCALE 

B2-B20+SCALE 

T»T0*SCAL£ 

FTR»0. 5/SQRT( AREA+AREAD) 
C"(B2^2,«T)>!'B]L*»2/2,-B2*(Bl-T)’^(Bl->-T)/2. 

C»C/AREA 

C 

C COMPUTE FACTORS FDR OIl/OV # DI2/DV * DJO/DV 
C 

0B1DW»B10*FTR 

DB2DV«B20*FTR 

DTDV=TO*FTR 

DIlD81»(B2 + 2.*T)’»' + 3/12.-B2<'*3/12. 

OIlDB2«(P2 + 2.*T)*<‘2*Bl/4.-B2**2*( Bl-T) / 4. 
DIlDT-(B2+2.*T)«*2*Bl/2.+32**3/12 , 

DC0B1«(B1*( B2+2.*T)-B2*S1-2.*C*T) /AREA 
DC0B2«(Bl^*2/2.-(Bl-T)*( Bl+T) /2.-C^r )/AREA 
DCDT»(B1**2 + B2«T-C«(2.’«'B1 + B2) )/AREA 
DI2D81a{T>*'81**2/2. ) + (2.<‘T*( B1/2.-C )**2)^ 

1 (4.*T*Bl«(Bl/2.-C)>!>( .5-DCDBl) ) + 

2 (2.*B21«T*(C“T/2.)*DCDB1) 
DI2D82*(-4.*T«B1«(B1/2.“C)’)'DCDB2 ) + ( r* + 3/l 2 . ) + 

1 (T»(C~T/2. )-«'*2) + (2.'«‘B2*T* (C-T/2. )*DCDB2) 

0I20Ta(Bl*>t'3/(S. ) + (2.>!<Bl>MBl/2.-C )>i'«2)- 

1 (^».>^T«Bl«(Bl/2,-C)*DCDT) + { R2«r**2/A.) + 

2 (B2t‘(C-T/2. )^«2} + (2.^B2*T*(C-r/2,)»(DCDT~.5) ) 
DJ0DB1=2 .’!'T*’!'3/3, 

DJODB2=T>f^3/3. 

0J0DT«(2.’4‘B1 + B2)*T’«'*2 
C 

C COMPUTE OIl/DV , DI2/DV , DJO/DV 
C 

0I1DV«DI1D81«DB1DV+DI1DB2«082DV+DI1DT*DTDV 
OI2DV»DI2DBl^r)RlDV + DI2DB2*D0 2DV+DI2DT*DTDV 
OJODV»OJODB1»'D81DV+DJODB2<'D32DV+OJODT*DTDV 
C 

C CREATE RUNSTREAM TO FIND DK/DV 
C 

WRITE(10,106)0APW,DUDV 
WRITE( 10# 107)DI2DV#DJ0DV 

106 FORMAT { 8X#*S l*SiiM( F13.5# 1X#«A1#»E 13,5# lX#tA2)*) 

107 fORMAT(8X#«S2 = SUM(<‘#E13.5# 1X# + A3#*E13.5# 1X#*AA) + ) 
WRITE{10#108) NDF#NDVJIM 

108 FORMAT(« DKDV SPAR ♦# 12# 14# ♦»SU<'U S 1 # S2 ) ♦ ) 

C 

C CREATE RUNSTREAM TO FIND DM/DV 


C 

WRITE(10#206 )DADV#DI1DV 
WRITE ( 10#207)DI2DV#DJ0DV 

206 FORMAT( 8X#*T1»SUM(«, E13.5# 1X#*B1# ♦E13.5#1X»*B2)^) 

207 F0RMAT(8X#*T2«SUM(*#E13,5# lX#*B3#tEl3.5#lX#+B4)*) 
WRITE ( 10#208 ) NDVJIM 

208 FORMATS DMDV DIAG 0*# I3#* = SUM('T1# T2 )♦ ) 
WRITE(10#209) 

209 FQRMAT<*CXQT DCU*/+ TOC !■«■) 

RETURN 

END 
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APPENDIX I.- LISTING OF BMSTRS 


SUBROUTINE BMSTRS USW? KC NT # J CNT» 1 B E AM ) 

COMMON KOREi>KEVEN?Z(lJ 

DIMENSION DY1DV(A),DY2DV(A) ,Y(A,2),S{8)>F3(2)»XM1(2)#XM2(2) 
DIMENSION X(20) 

NAMELIST/LINKF/NDV»X 
JMl - JCNT-1 
C 

C STORE FORCES AND MOMENTS IF ISTRS EQ 1 

C STORE DERIVATIVES OF FORCES AND MOMENTS IF ISTRS EQ 2 

C 

F3(l)»-Z( JMl+13) 

XMKl )— Z( JMl + lA) 

XM2(1)»-Z( JMl+15) 

F3{2)»Z< JMl+19) 

XMK2 )»Z( JMl+20) 

XM2(2)»Z( JMl+21) 

IF(KCNT.NE.l) GO TO 31 
IF(ISW.NE.I) GO TO 31 
C 

C STORE AREAf MOMENTS OF INERTIA, AND Y VALUES IF 
C FIRST TIME THROUGH 

C 

XIl»Z(JMl+25) 

XI2»Z( JMl+27) 

TCNT-38 
DO 3 K » 1,A 
DO 2 L » 1,2 
ICNT*ICNT+1 
Y(K,L)»Z( JMl+ICNT) 

2 CONTINUE 

3 CONTINUE 
C 

C FIND DERIVATIVE OF Y VALUES IF NEEDED 

c; 

c 
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C READ IN CONSTANT NUMBERS FROM UNIT 30 
C 

REAO(30»10) B10»B20>T0 
10 F0RMAT(3F10.3) 

READ IN DESIGN VARIABLES FROM UNIT 31 

READ(31#LINKF) 

CALCULATE DY/DV 

A0*(2. + B10-*-B20)«T0 
A*A0/X(2) 

$CALE»SQRT(A/AO) 

81»B10>>SCALE 
B2«B20»SCALE 
T»T0*SCALE 
FTRsO. 5/SQRT(AtA0) 

C»{B2 + 2, + T)*Bl>!<^2/2,-B2*(81-T)t(81+T)/2. 

C»C/A 

DB1DV«S10*FTR 
DB2DV»B20’»'FTR 
DTOV-TO+FTR 

0CDBl»<Bl*(8 2 + 2.*T)-B2'‘'Bl-2. *C*T )/A 
0CDB2»(B1^*2/2.-(B1-T)*(BH-T)/2.-C*T)/A 
0CDT=(B1**2 + 82*T-Ct(2 .♦Bl + 82) )/A 
DCDV*( DCDB1’«‘DB1DV)-MDCDB2*‘D82DV) + (DCDT^DTDV) 
DY1DV(1)»-DB1DV+DCDV 
0Y1DV(2)*DCDV 
DY1DV(3) «DCDV 
DY1DV(A)«DY1DV(1) 

DY2DV(1)=( .5*DB2DV) +DTDV 
DY2DV(2)»0Y2DV(1) 

DY2DV(3)«-0Y2DV{1) 

DY2DV(A)*DY2DV(3) 

31 ICNT»0 

DO 200 II = IfZ 
DO 100 I « l#4 
ICNT=ICNT+1 

COMPUTE STRESSES OR STRESS DERIVATIVES 

$(ICNT)=(F3(II)/A)+( (XMl(in/XIl)^Y(I,2) )-((XM2(II)/XI2) 
1 ♦Y(I,1) ) 

C 

C IF AREA OF BEAM IS A CONTRIBUTING FACTOR TO STRESS 
C DERIVATIVES AND ISTRS EQ 1 THEN CALCULATE THE FACTOR 
C AND STORE ON UNIT 16 
C 
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IF(ISW.NE.l) GG TO 95 

DFAC={(Xf^l(II)/XIl)<'GY2DVn))-((XM2(II)/;(I2)’«'DYlDV{I)) 
WRITE(16) DFAC 
GO TO 100 
C 

C IF AREA OF BEAM IS A CONTRIBUTING FACTOR TO STRESS 

C DERIVATIVES AND ISTRS £Q 2 THEN READ FACTOR OFF UNIT 16 

C AND ADD IT TO STRESS DERIVATIVE 
C 

95 IF( I8EAM.EQ.0) GO TO lOC 
READ<16) DFAC 
S( ICNT)«S( ICNT)+DFAC 
100 CONTINUE 
200 CONTINUE 
C 

C WRITE STRESSES ON UNIT 15 
C 

WRITE(15) S 

RETURN 

END 


9 
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